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Abstract 



The electroweak properties of nucleons and hyperons are calculated in a relativistic con- 
stituent quark model. The baryons are treated as three quark bound states, and the 
diagrams of perturbation theory are considered on the light front. The electroweak prop- 
erties of the baryons are of nonperturbative nature and can be represented by one-loop 
diagrams. We consider different extensions of the simplest model: 

• Quark form factors. 

• Configuration mixing of the wave function. 

• Asymmetric wave function. 

• Wave function different from the one of a harmonic oscillator valid up to energies 
of more than 30 GeV^. 

A comprehensive study of various baryonic properties is given: 

• Elastic form factors of the nucleon. 

• Magnetic moments of the baryon octet. 

• Semileptonic weak form factors. 

This analysis also gives the Kobayashi-Maskawa matrix element Vus and a sound sym- 
metry breaking scheme for the Cabibbo theory (see Sec. |4.4| ). 

A consistent physical picture appears in this work. The nucleon consists of an unmixed, 
symmetric three quark state, the wave function of the hyperons is however asymmetric 
with a spin-isospin-0 diquark. Only for the strangeness-changing weak decay do we need 
nontrivial form factors. 
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Chapter 1 
Introduction 



The purpose of this thesis is to present the results of comprehensive calculations of elec- 
tromagnetic and weak form factors of the baryon octet in a relativistic constituent quark 
model. Ever since the proposal of the quark model in the early sixties by Gell-Mann 
the modeling of the hadrons has been a very active area of theoretical research. There 
are many interesting current questions, which can only be answered by quark model cal- 
culations. New high-statistics data from hyperon beta decay raise questions about SU(3) 
breaking p, ^ and the related controversy about the meson and baryon derived values 
for the Kobayashi-Maskawa matrix element Vus has not been solved yet [^. This value 
is important for studying the unitarity of the Kobayashi-Maskawa matrix. There are 
also many recent papers on the magnetic moments of the baryon octet, which are a hard 
testing ground for various hadronic models. 

In this thesis we consider a relativistic constituent quark model on the light front 
that was first formulated by Terent'ev and Berestetskii 0, 0. It has been applied to 
various hadronic processes by Aznauryan et al. 0, ||, 0. Recently, new studies have been 
carried out by Jaus in the meson sector |I0| , [Tl|, |12[ and by Chung and Coester on the 



electromagnetic form factors of the nucleons |13|. Dziembowski et al. 1 14, 15, 16 1 and 
Weber et al. ||T^, treat the Melosh transformation of the quark spin in a "weak-binding 
approximation" of questionable validity. 

Nonrelativistic constituent quark models are successful in describing the mass spec- 



trum of baryons (for a review, see Refs. [|T9], |2y, ^). The dominant effects of the gluonic 



degrees of freedom are absorbed into constituent quark masses and into an effective con- 
fining potential. In addition, effective dynamics inspired by QCD has been considered 
in Refs. [p2| , |23| ; but this is inconsistent for light quarks. A review in Ref. p4| gives an 
estimate of the r.m.s. radius of the nucleon from 0.8 fm (derived from the charge radius) 
to 0.4 fm (obtained from hyperon decays). Considering the uncertainty principle these 
values of the r.m.s. radius imply a quark momentum in the range 250-500 MeV, which 
has to be compared with the light constituent quark mass in the range 210-360 MeV. 
The use of nonrelativistic quantum mechanics is therefore inconsistent, even for static 
properties of the hadrons, because the relativistic corrections are of the order (p^)/m^ 
and must be expected to be large. 

In a relativistic theory the Poincare invariance has to be respected; this means, on the 
quantum level, the fulfillment of the commutation relations between the generators of the 



1 



2 



CHAPTER 1. INTRODUCTION 



Poincare group. Dirac |2^ has given a general formulation of methods to simultaneously 
satisfy the requirements of special relativity and Hamiltonian quantum mechanics. An 
extension of the Dirac classes of dynamics can be found in Ref. . The light front scheme 



is in particular distinguished from the other Dirac classes. Among the ten generators of 
the Poincare group, there are in the light front approach seven (the maximal number) 
generators of kinematical character, and only the remaining three generators contain 
interaction, which is the minimal possible number. The light front dynamics is therefore 
the most economical scheme for dealing with a relativistic system. If we introduce the 
light front variables = ± p^, the Einstein mass relation p^p'^ = m? is linear in p~ 
and linear in p"*", in contrast to the quadratic form in p° and p in the usual dynamical 
scheme. A consequence is a single solution of the mass shell relation in terms of p~, in 
contrast to two solutions for p^: 



p ={pj_ + Tn)/p , p = ±Yp2 + . (1-1) 

The quadratic relation of p~ and p± = {p^iP^) in Eq. (|1 . 1|) resembles the nonrelativistic 
scheme [0, and the variable p"*" plays the role of "mass" in this nonrelativistic analogy. 
It is therefore a good idea to introduce relative variables like the Jacobi momenta when 
dealing with several particles. As in the nonrelativistic scheme such variables allow us to 
decouple the center of mass motion from the internal dynamics. Hence we do not have 
the problems with the center of mass motion which occur in the bag model. The light 
front scheme shows another attractive feature that it has in common with the infinite 
momentum technique [^]. In terms of the old fashioned (Heitler type, time ordered, pre- 
Feynman) perturbation theory, the diagrams with quarks created out of or annihilated 
into the vacuum do not contribute. The usual qqq quark structure is therefore conserved 
as in the nonrelativistic theory. It is, however, harder to get the hadron states to be 
eigenfunctions of the spin operator p9| . 



In this thesis we describe the baryon as a three quark bound state with relativistic Fad- 
deev equations and a Bethe-Salpeter interaction kernel. Using the quasiparticle method 
these equations can be reduced to a relativistic Schrodinger equation with an effective 
potential. We could either start with an Ansatz for the effective potential and derive the 
wave function or else start with the wave function itself. For simplicity we take the latter 
approach and work with an Ansatz for the wave function. In addition to the standard 
Gauss shaped wave function we also investigate a Lorentz shaped one that fits the nu- 
cleon form factors up to more than 30 GeV^. In addition to the minimal quark model 
(MQM), which uses the smallest number of parameters, we also consider some extensions 
to the MQM: structure of the constituent quark, several different radial wave functions, 
and asymmetry in the wave function. For the MQM we have five parameters: the quark 
masses m„ = = rriu/d and m^, and the wave function parameters j3 for the nucleons, 
the Ss and the Ss respectively, which essentially determine the confinement scale. For 
the asymmetric model we have three additional parameters: instead of (3 we have f3q and 
(3q. In the configuration mixing case there are some more parameters that determine the 
mixing and the confinement scale of the various radial wave functions. For the structure 
of the constituent quarks there are seven additional parameters, i.e. the anomalous quark 
moments /2m, /2d, /2s, and the vector and axial weak form factors fiud, fius, diud, and 
Qius at zero momentum transfer. Our goal is to keep the number of parameters as small 
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as possible in order to have a model with predictive power. On the other hand the MQM 
is not able to fit all electroweak properties of the baryons. The parametric dependence 
on the data is highly nonlinear, so it is not obvious that even 15 parameters may fit 
two experimental data. Our analysis gives a physical picture of a baryon, which is an 
asymmetric three quark state with quark form factors in the weak sector. 
This thesis is organized as follows: 

Chapter ^ In this Chapter we present the light front formalism for bound states of three 
quarks. We start with the light front variables, elaborate on the quasipotential 
reduction of the Bethe-Salpeter equation, and calculate the one loop diagram. We 
also study the wave function and compare it with the one given by Chernyak and 
Zhitnisky [0, and Lepage and Brodsky 



Chapter ^ In Sec. 3J- we give the details of our calculation of the electromagnetic form 
factors. We fit the mass rriu/d and the scale parameter (3 to the data of the magnetic 
moments of the proton and neutron, and the weak axial form factor gi{n p). We 
not only examine the Gauss shaped wave function, but also a Lorentz shaped one. 
With the latter it is possible to fit the data in the whole experimentally accessible 
energy region up to more than 30 GeV^. In Sec. |3.2| we derive the explicit formulae 



for the magnetic moments. The MQM is not able to give a reasonable fit, even with 
non-zero quark anomalous moments. Only an asymmetric wave function can fit all 
magnetic moments of the baryon octet as is shown in Chapter ^ The mass rus and 
the range parameters Pt./a and P~ can be fixed in Chapter |^. 

Chapter ^ Semileptonic beta decay of hyperons is an active area of current research, 
both experimentally with high-statistic data and theoretically with the problem of 
the Kobayashi-Maskawa matrix element V^s 0, 0]- In addition to the calculation 
of the weak form factors and their derivatives we present a Cabibbo fit that gives 
almost the same value for Vus as the one recently published [Q. 

Chapter ^ The smallest extension to the MQM is to use an asymmetric wave function, 
in which the scale parameter /? is replaced by two scales Pq and Pq corresponding to 
the quark-diquark binding and the quark-quark binding within the diquark. In the 
limit of two particles this extended MQM reduces to a MQM of the meson, which 
is successfully used in Refs. |1^ . 



Chapter ^ summarizes our investigation. We discuss the physical picture derived from 
our analysis and draw some conclusions. 

Appendix ^ describes in detail the symbolic and numerical methods used in this thesis. 

Appendix ?? contains all formulae for the spin matrix elements to all orders of the 
momentum transfer. 



The index at the end of this thesis should give quick access to special topics. 
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Chapter 2 

Light-front formalism for baryons 



2.1 Relativist ic three-body equation 

To specify the dynamics of a many-particle system one has to express the ten generators of 
the Poincare group and M^^ in terms of dynamical variables. The kinematic subgroup 
is the set of generators that are independent of the interaction. There are five ways to 
choose these subgroups [Q. Usually a physical state is defined at fixed xq, and the 
corresponding hypersurface is left invariant under the kinematic subgroup. 

We shall use the light-front formalism which is specified by the invariant hypersurface 
= + = constant. The following notation is used: The four- vector is given by 
X = {x^,x^,x±), where x^ = x^ i: x^ and x± = {x^yx"^). Light-front vectors are denoted 
by an arrow x = {x~^, x±), and they are covariant under kinematic Lorentz transformations 
The three momenta pi of the quarks can be transformed to the total and relative 



momenta to facilitate the separation of the center of mass motion p3 . 



P = P1+P2+P3, ^ = + , + , V= 

Pi + P2 



q± = (1 - OPi± - ^P2± , Q± = (1 - V)iPl± + P2±) - VP3± ■ 



(2.1) 



Note that the four-vectors are not conserved, i.e. Pi + P2 + Ps 7^ P- In the light-front 
dynamics the Hamiltonian takes the form 

H = . (2.2) 

2P+ ' ^ ^ 

where M is the mass operator with the interaction term W 



M = M + W 

= Oi ^M + J^L. (2.3) 

ri[l — rj) rj 1 — rj 
2 22 
Ml = + + ^ 

e(i-0 e 1-e 
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with rrii being the masses of the constituent quarks. To get a clearer picture of M we 
transform to gs and Q3 by 

Ei-\- q3 E12 + Q3 



El + E2 E12 + E3 

(2.4) 

i?i/2 = (q^ + m?/,)V2, E3 = (Q2 + m^)V2, = (Q^ + M3Y/^ 

where q = (gi, q2, gs), and Q = (Qi, <525 Qs)- The expression for the mass operator is now 
simply 

M = Ei2 + E3, Ms = Ei + E2. (2.5) 

We shall assume only two-particle forces interacting in a ladder-type pattern so that the 
dynamics of the three-body system is governed by the Bethe-Salpeter (BS) interaction 
kernel for the two body system and the relativistic Faddeev equations. 

Using the Faddeev decomposition for the vertex function F = F*^^^ + F*^^) + F'^^^ we can 
write down a BS equation for the various components in operator notation 

_ T(^)G2G3(r(') + F(3)) (2.6) 

with 

G,=:^,-m,, T« = (1 - VGaGa)-^^ , (2.7) 

and similarly for F*^^-* and F^^-*. V is the one gluon exchange kernel between two quarks, 
and T is already the ladder sum to all orders. It is useful to consider the second iteration 
of the vertex equation, which is given by: 

r = t/GiGsGgF , (2.8) 

where T = (F^^), F^^), F^^)) and U is the matrix 

^fr«G,rW for J with A;^^, J , . . 

\ T^'\GkT^'^ + GiT^''^) for i = j with k^l^i . ^ ' 



The four-dimensional Eq. ( ^.81 ) can be reduced to a three-dimensional equation 

r = WgsT , W = {1- URsy^U (2.10) 

by writing G1G2G3 = g3 + R3 where has only three-particle singularities. We choose a 
g^ which puts the quarks on their mass shells: 

g^ = {2mf [ ds-^ n 5+{p1 - m^){^i + m,) , (2.11) 

where P is the total momentum of the bound state, s = {P1+P2+P3Y and pi are restricted 
by pI > 0. We get 

,3 = {2mrS{pl - ml)S{pl - ml)Qm{l - O^ivMl - ^) (2-12) 



2. 1 . RELATIVISTIC THREE-BODY EQ UATION 
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with the spin projection operator 

A+fe)=E«fe'^)^fe'^)- (2-13) 

A 

Writing 

1 

93 



r« = Lppp Y^^u{vx\x)u{v2\2)u{v3\3) . (2.14) 



M^M'jMM' 



1/2 



we are led to the integral equation 

f»(q,Q,Ai,A2,A3) = E /ciVrfVW^'^(q,q,Q,Q',Ai,A;,A2,A'2,A3,A'3) 

x^3(q',Q')f(^)(q',Q',A;,A'2,A'3). (2.15) 

We can write this equation in terms of the wave function \E'. The Faddeev decomposition 
is \1/ = vl>(i) + _)_ v|/(3)^ the relation to the vertex function is = ^sr*^*^ and writing 

^ = (^(1)^^(2)^ ^(3)) 

(M^-M^)* = iy* (2.16) 

with Mb being the mass of the baryon. If we put W = MW + WM + W"^ we see that 
the wave function is an eigenfunction of the mass operator M^, given in Eq. ( |2.3| ): 

M^^ = Mlm (2.17) 



which is equivalent to the equation usually used in constituent quark models |[T9 



{E12 + E^ + W)^ = Mb^ . (2.18) 

This last equation is the starting point for an explicit calculation of the wave function, 
which has been done for the meson sector 
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2.2 Current matrix element 

We would like to calculate the current matrix element Af^ = {B'\q'~f'^q\B) corresponding 
to the Feynman diagram of Fig. |2.1| (a). 





+ 



+ 



(d) 



(e) 



Figure 2.1: Feynman diagrams for the elastic form factor of baryons. Only the three 
quark core of the baryon is considered. 



We restrict ourselves to the three quark core of the baryon. The diagrams of Fig. 
(c+d), and (e) can be absorbed into the wave function and quark form factors, respectively. 
We are left with the diagram Fig. ^]T] (b). 




+ 




+ 




+ 



(a) 



(b) 



(c) 



(d) 



Figure 2.2: (a) Feynman diagram for the elastic form factor, (b)-(d) Time x+-ordered 
diagrams corresponding to (a). Pointed lines represent instantaneous quark propagators. 



In Fig p.2| , the Feynman diagram (a) is equivalent to the sum of diagrams (b) - (d) 
in the old-fashioned perturbation theory. If we consider x^-ordering and put = 0, 
diagram (c) drops out because of conservation of +-momentum. The x^-instantaneous 
propagator in diagram (d) is proportional to 7^, and gives no contribution for = 



2.2. CURRENT MATRIX ELEMENT 
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since (7''')^ = 0. We are left with diagram (b) which can be expressed with the help of 
the vertex function F [Nc is the number of colors): 



'2 



indices V^l ~ "'1 l/^l " "^1 + «ej 

ki 



xFt ^'(P2A2)n°(p2A2)^ (p3A3)nnP3A3) 

(pi - mi + ze) (pi - mi + ^e) + P^^^tions , 
Pi = \p + k + \l, p2 = \p-k + ]4, p, = ^P-r, (2.19) 
pi' = pi-K, K = P-P' . 

On the light front, we have exact correspondence with the choice of the Greens function 
in Eq. ( |2.12| ). If the vertex function F is assumed to be independent of the components 
k~ and /~, we can calculate M"*" by contour methods in the k~ and /~ planes. M"*" is 
given by the residua of the two noninteracting quark poles. Replacing vertex functions 
by wave functions we get 



x(Oi + 02 + 03)^(q,Q,Ai,A2,A3 
1 



Oi = -^m(pi'A;)7MpiAi) , (2.20) 



O2 = —-^——u{p2'X2h^u{p2X2), 
^3 = (1 _ ^^ ^(P3'A3)7+m(p3A3) . 

The OjS correspond to the zth diagram in Fig. |2.3| . For Oi the primed variables are 

q^ = q±-il-OK±, Q'^ = Q^-il-r])K^, (2.21) 

for O2 

q'^ = q^+^K^, Q'^ = Q±-{1-V)K^, (2.22) 

and for O3 

q'± = q±, Q'^=Q±+VK±. (2.23) 
For pointlike quarks the matrix element of the current is 

u{p^'X[h+u{p,Xi) = 2^7]P+6x,y^, 

uip2'X',)j^u{p2X2) = 2(l-Or7P+5A,A^ , (2.24) 
u{P3'KhMP3h) = 2(l-77)P+5,3,, . 

After factoring out color the wave function is totally symmetric and we have J2i (Oi) = 
3 (Oj) for any j. Since primed variables take a simple form for O3 we choose 3 (O3). We 
arrive at the form 

1 /2 

= 2P+^ / d'qd'Q (flf^) E ^^(q'. Q'. A'3)35A3A^^(q, Q, A3) . (2.25) 
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Figure 2.3: Feynman diagrams that represent the transition of the baryon state with four- 
momentum P to the baryon state with four-momentum P'. K = P — P'. The photon or 
the W boson is coupled either to the first, second or third quark hne, corresponding to 
the diagrams (1),(2) and (3), respectively. 



2.3 Wave function 

In light front variables one can separate the center of mass motion from the internal 
motion. The wave function is therefore a function of the relative momenta q and Q. 
The product = ^x^P^ with $ = flavor, x = spin, and (j) = momentum distribution, is a 
symmetric function. This is consistent with Fermi statistics since the color wave function 
is totally antisymmetric. 

The angular momentum j can be expressed as a sum of orbital and spin contributions 



j = «Vp X p + , 



^2.26) 



where TZm is a Melosh rotation acting on the quark spins s^, which has the matrix repre- 
sentation (for two particles) 



{\'\nM{^,q±,m,M)\\) 



m + — za ■ (n X q) 
m + ^M)2 + ql 



(2.27) 



A' A 



with n = (0, 0, 1). In previous works |T^, [T^, |1^, |T8[ this rotation has been approxi- 
mated by putting M = Mb- This corresponds to a weak-binding limit which cannot be 
justified for a bound state in QCD. In this limit our model has a close connection to many 



other relativistic quark models as shown by Koerner et al. [37 



The operator j commutes with the mass operator M; this is necessary and sufficient 
for Poincare-invariance of the bound state. 

In terms of the relative momenta the angular momentum takes the form 



j = iVq X Q + 7^M(r/,Q±,M3,M)jl2 + 7^Af(l-r/,-Qx,m3,M)s3 

ji2 = iVq X q-^7^A^(^,g±,^^l,M3)sl -^7^A^(1 -g±,m2,M3)s2 . 
We can drop the orbital contribution. 



(2.28) 



2.3. WAVE FUNCTION 
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7^l 



7^2 



7Z, 



ac - quQi 
cQr + aqR 



-aqi - cQl 
ac - qiQn 



ad + q^Qi aqi-dQi 
(P + g2 V dQa - aqR ad + qiQa 

h 



(2.29) 



-Q 



R 



with 



a 
c 
qR 
Qr 



mi + 

qi + iq2 , 
Qi + iQ2 , 



Ql 

b 



6 = + (1 -r7)M , 
, = m2 + (1 - 
qL = qi- m , 
Ql = Qi- iQ2 ■ 



(2.30) 



The momentum wave function (f) is normahzed according to 



d^qd^Q\4>[^ 



(2.31) 



and can be chosen as a function of M to fulfill the spherical and permutation symmetry. 
That is the same as to express it in terms of the off-shell energy S since 



The S-state orbital function 0(M) is approximated by either 



0(M) = A^exp 



M2 



or 



0(M) 



N' 



(M2 + /32)3-5 



(2.32) 



(2.33) 



which depend on two free parameters, the constituent quark mass and the confinement 
scale parameter j3. The first function is the conventional choice used in spectroscopy, but 
it has a too strong falloff for high K^. Both functions give nearly the same result for low 
values of K^, the second one performs obviously better for high K^. This independence 
of the wave function for low suggests that the static properties are mainly given by 
the flavor and spin part of the wave function. 

The total wave functions for the baryon octet are ^ 



P 



n 



A 



(uudx^^ + udux'^'^ + duux^^ 
v3 ^ 

^ (ddux^^ + dudx^^ + uddx^^ 
v3 ^ 

-1 
7^ 

^^{uusx'^ + usux'^ + suux'^)^ 



{uds — dus)x'^^ + {usd — dsu)x''^ + {sud — sdu)x 



p2 



;2.34) 



^The overall sign for E+,A and .n^ has to be changed in Ref. 



p. 46]. 
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with 



■^0 



-1 r 

-1 



{uds + dus)x^^ + {usd + dsu)x^'^ + {sud + sdu)x^^ (t> , 



71 



,A3 



\ 2 X ] 



^ {ssdx^^ + sdsx^^ + dssx^^ 



V3 



0, 



.A3 



xf 



^UTT + TiT -2Tn), 
^(2 iit - iti - Tii) 

^(Tit - itt) , 
^(Tii-in)- 



(2.35) 



The spin wave functions x^^ X^^ ^^e the appropriate permutations of x^^, and x^'^ 
and x''^ s-re the appropriate permutations of x''^- The spin- wave function of the ith quark 
is given by 

(2.36) 



T=7e,Q) and • 
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2.4 Extensions of the model 

As already mentioned in the introduction we are also considering some extensions to the 
minimal quark model (MQM), because the MQM is not able to fit experimental data for 
both the electromagnetic and weak sector with the same parameters. We give an overview 
of the different models in this Section. Questions concerning particular experimental data 
are discussed in the appropriate Chapters. 



2.4.1 Minimal quark model 

The MQM uses the wave function presented in Eq. ( p.33| ) with structureless quarks. It is 
called minimal because it uses the smallest number of parameters. These are the mass of 
the up and down quark m„ = rrid = rriu/d, the mass of the strange quark m^, and the wave 
function parameter j3 for the nucleons, the Es and the Ss; this (3 essentially determines the 
confinement scale. With these five parameters we either fit the electromagnetic properties 
(parameter set 4 of Table |3.1j on page ^) or the semileptonic weak decays (parameter 
set 6). The contrary statement in Ref. has to be questioned, since their numerical 
results for the magnetic moments of the baryon octet are wrong. 



2.4.2 Quark structure 

One extension to the MQM is the introduction of nontrivial quark form factors. This 
would give us three new parameters in the electromagnetic sector (quark anomalous mag- 
netic moments), and four new parameters in the weak sector (axial and vector quark 
form factors at zero momentum transfer). Chung and Coester investigate the nucleon 
sector of the same model, and favor quark anomalous magnetic moments and a modified 
axial coupling for quarks. There are two reasons why we think that the quarks should be 
structureless in the electromagnetic sector: 

1. Whatever the nature of this form factor is, whether it be due to a composite model 
or to radiative corrections, one expects the quark form factors to fall off for large 

in a different way from that used in Ref. |TB[. 

2. There exists no parameter set of the anomalous magnetic moments that can improve 
the magnetic moments of the baryon octet (see Sec. 



2.4.3 Configuration mixing 

Another extension of the MQM is the admixture of different radial wave functions. The 
configuration mixing suggested by spectroscopy reads: 

iBaryon) = A [56,0+] [8 x 2] + B [56,0+]*[8 x 2] + C [70,0+] [8 x 2], (2.37) 



in the notation [SU(6),LP][SU(3)flavourXSU(2)spin], where A^ + B^ + C^ = 1 , L denotes the 
angular momentum, and p is the parity of the nucleon. The values for A, B, C are listed 
in Table Ol for different references. 
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Table 2.1: Parameters for the configuration mixing of the baryon octet given in Eq. ( |2.37| ) 
for two different references. 







A 


B 


C 


Ref. 


40 


0.93 


-0.29 


-0.23 


Ref. 


m 


0.90 


-0.34 


-0.27 



In principle no additional parameters are required, if we use the same scale parameter 
P for every wave function, but in practice, it is convenient to choose different f3s or 
admixture parameters. The wave functions in Eq. (|2.37|) are as follows: 

[56,0+][8x2] = -1= (x^<i>P + x'<i>'' 



[56,0+]*[8x2] = 1^ (x^c^P + x''^') € , (2.38) 



V2 

[70,0+][8x2] = Ux'^' + x'^'')<l>p+Ux'^'-x'^ 



The spin functions x'' ^^'i X'^ ^-I's the same as in Eq. ( |2.35D , the flavor wave functions 
and $^ correspond to x'^ X'^ with spin up and down exchanged with the appropriate 
flavors of the baryon. For the momentum wave functions (ps, <P*si 4'p, and (px we first define 

M, = ^^^^, Xi = pt/P+, h^=Pi±-XiP^. (2.39) 

With these functions it is easier to build wave functions with special symmetries. Note 
that M = Ml + M2 + M3, the combination Mi + M2 — 2M3 is symmetric in the particles 
(12), and Mi — M2 is antisymmetric in the same particles. We therefore write: 

CPs = N,e-^'/'^\ 

0: = iV:(MV/3'-c)</>,, 

0A = iVA(Mi + M2 - 2M3)0. , 

0, = Np{Mi-M2)(Ps. (2.40) 
The constant c in (p* is evaluated from the orthogonality of (ps and (p*: 

and the constants A^^, A^*, Nx, and Ap are given by the normalization in Eq. ( p.31| ). 
These wave functions in Eq. ( |2.40| ) go over into the nonrelativistic ones ||T, 3U| in the 
limit, where the masses rrii go to zero and the ^ and r] to their nonrelativistic values. 

Unfortunately, the mixing configuration does not improve the fit, it is even worse for 
the crucial ratio in Eq. ( |4.27] ). A rough estimate gives 

gi//i(A ^ pe-^J ^ _3 A 8 ^ ^ ^ _ 

gi/fi{T.~ ne-z/e) V 3 / 
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to be compared with the MQM value —3, and the experimental data —2.11 ±0.15. Other 
values like the ratio fi{p)/fi{n) also get worse. We therefore do not consider this extension 
any further. 

2.4.4 Asymmetric wave function 

The extension to the MQM with a two quark clustering in the light front wave function 
is minimal in the sense that there is no difference between the extension and MQM in 
the meson sector. Instead of one scale parameter /3, we have two of them, f]q for the scale 
between the two spin-isospin-zero quarks and Pq for the scale between the third quark and 
the diquark. We devote the entire Chapter ^to the asymmetric wave function, because it 
improves the fit dramatically in many details and it provides a comprehensive fit of both 
the electromagnetic and weak sectors. 
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2.5 Discussion of the wave function 

An analysis of the baryon spectrum based on Eq. ( p.l8|) could in principle determine 



the wave function, but we restrict ourselves to the two approximations in Eq. (|2.33|) . It 
is therefore important to compare our Ansatz with the wave function of other authors. 
Usually the transverse momenta are integrated out up to a scale fi, and the wave function 
is expressed in light front fractions Xi = pf /P^, written in our variables as 

xi = ^r] , 

X2 = ^7(1-0, (2.43) 

X3 = l-T] . 

The valence quark distribution amplitude ^(xi, 0:2, X3, /x^) is 

/•kil<A'^ r\Ql\<t^^ 

^{xi,fi^)=J J ^{x,,q^,Q^)d^q^(fQ^ . (2.44) 

This amplitude is well known in two limits, which are unfortunately not interesting. In 
the static, symmetric SU(6) quark model, the variables Xj take on only discrete values: 

0nr(x.) = 5 (xi - ^) 5 (x2 - ^) (^3 - ^) = 5 - ^) 5 (r^ - ^) (2.45) 



and the asymptotic amplitude 0as for large K'^ is known as 

4>-as{xi) = (l){xi,fi^ ^ 00) = 120x1X2X3 , (2.46) 
with a normalization such that 

J 0(xj, /i^)5(^ Xj — I)dxidx2dx3 = 1 . (2.47) 

Notice that this is not the same normalization as the one used in Eq. ( |2.31| ). Unfortunately 
the knowledge of both forms is not very useful since they contradict experimental data 

111- 

Using the QCD sum rule technique, Chernyak and Zhitnitsky |30| suggest the quark 



distribution amplitudes for the proton as follows: 

0cz(a;i,/i^ 1 GeV) = (2.48) 
120x1X2X3 [ll.35(x? + xl) + 8.82x^ - 1.68x3 - 2.94 - 6.72(x^ - xj) . 
For the Gauss shaped wave function (pQ and the Lorentz shaped one 0l in Eq- (p.33|) we 



can write 



(MV/32 + I)" 
/ d'q^d^Q^ = 

-00 J— CO 



NP',„HI - ,)(1 - ?) exp - - j . (2.49) 



/OO roo 
/ 0L d^^d^Qj 



I —00 J —00 

iV'/52'^(l - r7)"~ir/"(l - ^)'^-i^"-i 

K(i - 0(1 - ^) + ^ie(i + - + /5^^e(i - 0(1 - V)] 



n-2 



2.5. DISCUSSION OF THE WAVE FUNCTION 



17 



Letting the quark masses m, go to zero the amphtudes both converge to the asymptotic 
form (|2^) : 

0g(x.) > iV/5V(l-^) (1-0=^/3^X2X3, 

M^^) ^ N'P^^r]\l-T])il-0=N'P^x,X2X,. (2.50) 

The differences between these various wave functions are best seen in a plot. In Fig. p.4| 
we show 0as (a), 0cz (b), 0g (c+d), and 0l (e+f). The plots (c) and (e) are the symmetric 
wave functions (parameter set 6), (d) and (f) are the asymmetric ones (parameter set 8 for 
hyperons). The broad, unstructured distribution in the asymptotic limit gets sharper and 
more structured for the phenomenological amplitudes. The wave function in (c) usually 
used in quark models is close to the asymptotic function (a). 

The important difference between (pc and 0l is their large momentum behavior. For 
\K'^\ — >• oo the wave functions behave as: 

(2.51) 

The exponential falloff for (pc becomes too strong at a momentum scale of about 2 GeV^ 
(see Fig. 



5g 
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Figure 2.4: Different wave functions used by other authors and used in this work, (a) 
asymptotic quark distribution (b) amphtude derived by QCD sum-rule technique 
pO| ; Gauss shaped wave function with parameter set 6 (c) and set 8 for hyperons (d); 
Lorentz shaped wave function with parameter set 6 (e) and set 8 for hyperons (f). A 
trend to more structured, asymmetric wave functions can be seen. 



Chapter 3 

Electromagnetic properties 



3.1 Nucleon electromagnetic form factors 
3.1.1 Calculation of nucleon form factors 

The electromagnetic current matrix element can be written in terms of two form factors 
taking into account current and parity conservation: 



(AT, X'p'\J'^\N,Xp)=Uy{p') 



2M, 



N 



(3.1) 



with momentum transfer K = p' —p and = q'j'^q. For = the form factors Fi and 
F2 are respectively equal to the charge and the anomalous magnetic moment in units e 
and c/Mat, and the magnetic moment is = Fi{0) + ^2(0). The Sachs form factors are 
defined as 

Ge = F, + -—^F2, and Gm = Fi + F2, 



and the charge radii of the nucleons are 
,dFi{K'^ 



6- 



and 



E/M 



6 



dGE/M{K'^) 



dK^ 



K2=o ^E/AdiO) 

The form factors can be expressed in terms of the + component of the current: 

1 



(3.2) 



(3.3) 



K^F2{K^ 



2P 

2P+ 



iV,T 



(3.4) 



(iV,T |j+|iv,i) 



Therefore Eq. ( p.25| ) can be used to calculate the form factors. In addition, the Dirac 
quark current q'~f^q can be generalized to include the quark structure in the following way: 



(3.5) 



19 



20 



CHAPTER 3. ELECTROMAGNETIC PROPERTIES 



We get 



(2vr) 



d^qd'^Q 



E3E12M' 



1/2 



0t(M')0(M) 



3 r 



-2M. 



N 



(2vr)' 



(3.6) 



1/2 



0"f(M')0(M) 



X 



1=1 



with 



-1 

1 



i = (uud) for the proton and i = (ddu) for the neutron. The form 

factors Fiu{0) and Fid{0) are the charges, F2u{0) and -^2^(0) are the anomalous magnetic 
moments of the u and d quarks, respectively. Only the Fij-terms contribute for K'^ = 0, 
and <^Xt*IXt*) = 1- The matrix elements ^Xt*lxJ'*) ^-iid ^Xtl^slx^*) given in the next 
Subsection. The expressions for K'^ 7^ are quite long, we therefore give only Fi for the 



proton with vanishing quark anomalous moment and Fi^ 

Num 



'Id 



1=1 



(3.7) 



Num = [d^ + Q'l)[a^ + Ql){h'h + Q'^-Q^){^d^ + q'i) 

+(c2 + d^)ql{a'a\W + g^-Qx) + («« + ^QrQ±) 
x[266'(Ql-g^) + + {Q'±-Qa-?]} 

+cdql[Aaa'hh\Q'^-Q^) - 2{a^Q'l + a^Ql){hh' + Q'rQ±) 

For i^'^ = 0, Eq. ( ^.7[ ) reduces to 1 and we get -Fip(O) = 1, the charge of the proton in 
units of e. 

3.1.2 Results and conclusions 

Appendix ^ describes how the formulae were generated and integrated. The exponential 
function 0(M) in Eq. (|2.33| ) falls off too fast, it can only be valid for low K^. In general 
has to be just a function decreasing with M. We try 



(3.8) 



3. 1 . NUCLEON ELECTROMAGNETIC FORM FACTORS 



21 



and is chosen so that Eq. ( |2.31| ) is fulfilled. The wave function <^(M) does correspond 
to a confining potential in the sense that we need more energy to ionize the bound state 
than to produce a new quark pair. This guarantees that no free quark appears. Figures 
37^ and p.3| show that the low K"^ behavior is the same for both wave functions. 

We plot the parameter (3 against m^/^ for the different experimental data. In Fig. ^.1| 
the fits for the Gauss shaped wave function (a) and for the Lorentz shaped one (b) are 
given. We see that it is not possible to exactly fit the three data, since the three lines 
should ideally meet in one single point. Figure |3.1| shows for case (a) that the proton and 



0.8 



0.6 



0.4 



0.2 





\ l^(P) 


§1 


|x(n) 






(a) 







0.17 0.21 0.25 0.29 0.33 



1 




\ H(P) 




0.8 






gi 


0.6 








0.4 








0.2 


(b) 







0.21 



0.25 



0.29 



m, 



u/d 



m 



u/d 



Figure 3.1: The lines represent a set of parameters P and rriu/d, which reproduce respec- 
tively the experimental data for the magnetic moments of the proton and neutron, and 
for gi{n — >■ pe~z/). (a) Parameters for the Gauss shaped wave function; (b) parameters 
for the Lorentz shaped wave function. 



neutron magnetic moments alone tend to a large quark mass mu/^ — 0.33 GeV, whereas 
Qi of the neutron beta decay together with the neutron magnetic moment favors a small 
rriu/d — 0.23 GeV. We compromise and fit to the proton magnetic moment and to gi of 
the neutron decay , which yields 171^/^ ^ 0.267 GeV. The anomalous moments of the u 
and d quarks are fitted to F[{0) of the neutron. To analyze the results we have chosen two 



sets of parameters given in Table |3.1| on page ^ set 1 with quark anomalous magnetic 
moments and set 2 without them. The situation for case (b) in Fig p.l| is similar, with 
the exception that the proton and neutron magnetic moments do not favor a large mass. 
We use parameter set 3 for the Lorentz shaped wave function. 

Figures P]^ - E]^ show Fi and F2 for the proton and neutron for the various versions 
of the model. If we neglect the quark anomalous magnetic moments, only the form factor 
Fin (Fig W^ ) changes and (Fig \i.7\) gets shifted by a small amount. Both changes are 
welcome but only of minor importance. Fin is very small and therefore sensitive to any 
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Table 3.1: The parameters of the constituent quark modeL All numbers are given in units 
of GeV. Note that only set 3 is used for the Lorentz shaped wave function. 





rriu 


rrid 


rris 


Pn 








F2d 


F2s 


Set 1 


0.267 


0.267 




0.56 






-0.0069 


-0.028 




Set 2 


0.267 


0.267 




0.56 






0.0 


0.0 




Set 3 


0.263 


0.263 




0.607 






0.0 


0.0 




Set 4 


0.33 


0.33 


0.55 


0.16 


1.00 


1.08 


-0.0086 


-0.034 


0.077 


Set 5 


0.267 


0.267 


0.33 


0.56 


0.63 


0.70 


0.0 


0.0 


0.0 


Set 6 


0.267 


0.267 


0.40 


0.56 


0.60 


0.62 


0.0 


0.0 


0.0 


Set 7 


0.267 


0.267 


0.40 


0.56 


0.60 


0.62 


-0.0069 


-0.028 


0.056 



Table 3.2: The quantity -F'(O) for the nucleons with parameter sets 1 and 2 of Table [311 . 
The values for the Lorentz shaped wave function (set 3) are almost the same as the one 
from set 2. 



Form factor Set 1 [fm^] Set 2 [fm^] Expt. [fm 



F' 
F' 

^2p 

F' 
F' 

^2n 



0.0874 
0.179 
0.0027 

-0.186 



0.0924 
0.177 
0.012 

-0.170 



0.0966 
0.234 
0.0017 

-0.236 




Figure 3.2: The proton form factor Fip{K'^). Continuous line, Parameter set 3; dashed 
line, parameter set 2; dashed-dotted line, parameter set 1. The experimental points are 
taken from Ref. g^, |3 . 
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(N 




0.5 0.75 

K2 [GeV2] 



Figure 3.3: The proton form factor F2p{K^). Continuous line, Parameter set 3; dashed 
hne, parameter set 2; dashed-dotted hne, parameter set 1. The experimental data are 
taken from Ref. E|, ||, |5 



corrections. 

All data for form factors can be reproduced well up to 2 GeV^ for both wave functions 
in Eq. ( p.33 ). But only the Lorentz shaped wave function has a good high energy behavior 
up to more than 30 GeV^. In this region we already have QCD predictions for Gm ^ • 



Note that this is an extremely large K region, since other models can only fit the data 
up to 2.5 GeV^ ^ 



4^ and 6 GeV |T^. Our fit is even better than the well known dipole 



formula 



G 



Mp 



1 



-2 



Mv = OM GeV 



(3.9) 



One can approximate our form factors with the parameterization of Eq. ( [4.23| ) on 
page It is valid up to 3 GeV^ with a deviation of less than 5%. This justifies the use 
of that parameterization in Chapter ^. 

The values for magnetic moments derived with the parameter sets 1 and 2 have been 
collected in Table ^]3| and will be discussed further in Sec. |3]^. We do not consider the 
Lorentz shaped wave function (set 3) any further in this thesis, because the results for 
small momentum transfer or even static properties are almost the same for both types of 
wave functions. 
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-K2 [GeV2] 



Figure 3.4: The proton form factor Fip{K'^). Continuous line, Parameter set 3; dashed 
hne, parameter set 2; dashed-dotted hne, parameter set 1. The experimental data are 



taken from Ref. E|, 44, 45 




Figure 3.5: The proton form factor F2p{K^). Continuous line. Parameter set 3; dashed 
line, parameter set 2; dashed-dotted line, parameter set 1. The experimental data are 
taken from Ref. ra, M, H5 
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i 2 3 4 



-K2 [GeV2] 

Figure 3.6: The neutron form factor Continuous line, Parameter set 3; dashed 

hne, parameter set 2; dashed-dotted hne, parameter set 1. 




Figure 3.7: The neutron form factor F2n{K'^)- Continuous hne, Parameter set 3; dashed 
hne, parameter set 2; dashed-dotted hne, parameter set 1. The experimental data are 
taken from Ref. g|, ||, ig. 
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-K2 [GeV2] 

Figure 3.8: Proton form factor GMpiK"^) for both wave functions in Eq. |2.3tj| . Continuous 
line, Lorentz shaped wave function with parameter set 3; broken hne; Gauss shaped wave 
function with parameter set 2. The experimental data are taken from Ref. ||5T|| . 
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3.2 Magnetic moment of the baryon octet 
3.2.1 ^2(0) in the quark model 



According to Eq. (O) the magnetic moment of a baryon is fi = -Fi(O) + -F2(0). The form 
factor -Fi(O) is equal to the charge, and F2{0) to the anomalous magnetic moment k of 
the particle. We have [see Eqs. ( p.25| ) and ( |3.5| )] 



= 2P+^ / d^d^Q (flf^) ^^q, Q', A')OM/(q, Q, A) , (3.10) 
where O is given by: 

We get for the baryon octet 



i=l 



—1 

with O.s = ( Q ) , and with 



(3.11) 





P 


n 


S+ 






A 




■^0 


E^A 


a/b 


A/A 


A/A 


A/A 


A/A 


A/A 


p/p 


A/A 


A/A 


p/A 


i=l 


u 


d 


u 


d 


(u+d)/2 


(u+d)/2 


s 


s 


(d-u)/2 


i=2 


u 


d 


u 


d 


(u+d)/2 


(u+d)/2 


s 


s 


(d-u)/2 


i=3 


d 


u 


s 


s 


s 


s 


d 


u 





In order to get -Fi(O) and ^2(0) the matrix elements {xilxi) (Xtl^slXi) have to 
be calculated to order 1, and (xtIxj.) to order K±: 



(xflOalxf) = {xf\Os\xf) = 0, 
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Xf\0s\xf 



b^ + Q 



2 ' 



Xflxf) = Ik^{2A^ + 2A2-As) 



Xflxf) = \K^i2A, + 2A,-A2), (3.12) 



3 

Xflxf) = ^^^(2^2 + 2^3-^), 



Xflxf) = K^A 
(xflxf) = K^A, 

oil Al\ 1 



3 



2 



with 

Ai 



12 



2 ' 



^3 



+ Q 

^ - 2(ij^ ; £ 

^(^-2(1^) 

a? + Ql c2 + gi ■ 



Eq. (|3.10|) together with Eqs. ( p.ll| ) and ( p.l2| ) is the most general formula for ^2(0). 
Putting niu = rrid for the hyperons there is an equality A2 = A^ under the integral, which 
reduces the number of integrations. The simplified formulae read: 



~ 3^^ +3^^ +3 Um„ 2mJ 

I \ 2 (3) 2 (2) 1 / F2U 4:F2d\ 

'^^"^ = 3^-3^-3^-2;;^;^"' 

^ ' 3 ^ 3 ^ 3 2m„ ^ 3 2m, ^ ' 

^ ' 3 ^ 3 ^ 3 2m, ^ 3 2md ^ ' 

^ ^ 3 ^ 3 ^ 2m, ^ ' 

«:(--) = -Ui'^ --S^^ --^Z^i'^ + -^Z^^\ 
^ ' 3 - 3 " 3 2mrf " 3 2m, " ' 

«:(S°) = ^^4=^) - ^54^) - 1^Z4=^) + , (3.13) 

^ ' 3 " 3 ^ 3 2m„ - 3 2m, - ' ^ ' 
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k(S°A) 
with 

q{2) 



V3 



{S. 



(2) 
E 



A 



(2) 



5 



(3) 
B 



A 



(2) 



(3) 



-2Mb 
-2Mb 
-2Mb 
-2Mb 
2Mb 



{2ny 
{2ny 



(2vr) 



The masses rrii are set as follows: 
0(2) q(3) 7 

s'^\At\z§^ 

0(3) ^(3) 7(3) 
^^2) ^{2) 



J_ / _ F2d \ ^(2) 

v/3V2m„ 2mJ ^ 



rfV'Ql0l'(2Ai + A2)/3 
£qd^Q\(l)\\AA2 - Ai)/3 
d'qd^Q\(j)\''A2 , 



&2 + g2 



mi 
mi 
mi 
mi 
mi 



m2 
ms 
m2 
mg 
m2 



mg : 

m, 

m, 

m,, 

m,. 



m 



m 



m2 

ms = m^ , 
m2 = m , 
m3 = m , 



where m 



m„ 



m^. In k(S°A) the mass is equal to (Ma + Mj^)/2. 



(3.14) 



3.2.2 Results and conclusions 



Table |3.3| shows the magnetic moments for all baryons and the transition moment TP — >■ 
A7. For comparison we have chosen set 4 from Table 

These parameters are fitted to the magnetic moments and good agreement can be 
achieved even without anomalous quark magnetic moments. But we have 



(3.15) 



which cannot be reconciled with the weak beta decay, since the small wave function 
overlap would cause too large a suppression of the A^* = 1 transitions. Therefore, we 
reject the parameter set 4. In Chapter ^ it is shown that beta decay requires 



/?jv ^ /3s ^ /3e 



(3.16) 



So another set of parameters (set 6) will be considered in Chapter ^, but Table ^ shows 
that it leads to results which are in disagreement with experiment. Use of set 6 with 
non-zero quark anomalous magnetic moments (set 7) leads only to a slight improvement. 
For the quark form factor ^2^ we get 0.056 GeV by fitting /i(A). 

We did not used the Lorentz shaped wave function for the magnetic moments because 
this wave function is only important for high i^T^-values. 
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Table 3.3: Magnetic moments of the baryon octet and transition moment for S A7 in 
units of the nuclear magneton. 



Particle 


Experiment |52| 


Set 4 


Set 6 


Set 7 


P 


2.79 ±10-^ 


2.85 


2.78 


2.78 


n 


-1.91 ±10"^ 


-1.83 


-1.62 


-1.73 


S+ 


2.42 ± 0.05 


2.59 


3.23 


3.15 


s- 


-1.157 ± 0.025 


-1.30 


-1.36 


-1.56 


A 


-0.613 ± 0.004 


-0.48 


-0.72 


-0.58 




-1.250 ± 0.014 


-1.25 


-1.87 


-1.64 




-0.679 ± 0.031 


-0.99 


-0.96 


-0.71 




1.61 ± 0.08 


1.22 


1.74 


1.79 



We conclude that, within the symmetric wave function model, either the magnetic 
moments can be fitted and the weak decay parameters are poorly fitted or vice versa. 
The opposite statement in Ref. has to be questioned, because their numerical results 



for the magnetic moments are wrong. Our results agree with Ref. |5^ on this point. 
The inconsistency just described between the electromagnetic and the weak sector can be 
resolved by using asymmetric wave functions as described in Chapter ^. 



Chapter 4 

Hyperon semileptonic beta decay 



There is evidence for SU(3) symmetry breaking in semileptonic beta decay of hyperons 
0, Up to now one does not have any means to calculate such a behavior from first 
principles. Some symmetry breaking schemes have been proposed in the framework of 
bag models |^ and chiral perturbation theory |Q. Since they are different from each 
other we need some support from other models. 

Another uncertainty in the analysis of data is the dependence of the form factors. 
Although it is generally small, a change of My or by ±0.15 GeV in the case E~ nev 
(which has the largest -ft'max) causes a relative change of gi/ fi of ±2%. Ignoring the K"^ 
dependence altogether would shift gi/ fi by 17%. Our quark model provides a unique 
scheme for the calculation of these form factors. 



4.1 Hyperon semileptonic decay 

In the low energy limit the standard model for semileptonic weak decays reduces to an 
effective current- current interaction Hamiltonian 

H,nt = -^J^.L^ + \^.c., (4.1) 

where G ~ 10^^/Mp is the weak coupling constant, 

L^^ = V^e7^(l - 75)V^. + ^"m7^(1 - 75)^. (4.2) 

is the lepton current, and 

J,, = V^,- , V^ = VudUlt^d + VusU'j^s , Af, = VudU'j^'j^d + KsM7m75S , (4.3) 

is the hadronic current, and Vud, Vus are the elements of the Kobayashi-Maskawa mixing 
matrix. The r-lepton current cannot contribute since m,- is much too large. 
The matrix elements of the hadronic current between spin-^ states are 

{B',p'\V^^\B,p) = Vg,>u{p') 



Mi 



Mi 



u{p) , (4.4) 
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{B\p'\A^\B,p) = V,,,u{p'[ 



Mi 



Mi 



75m(p) , (4.5) 



where K = p—p' and Mj is the mass of the initial baryon. The quantities f\ and gi are the 
vector and axial-vector form factors, /2 and g2 are the weak magnetism and electric form 
factors and /a and 5^3 are the induced scalar and pseudoscalar form factors, respectively. 
T invariance implies real form factors. We do not calculate /s and since we put K'^ = 
and their dependence on the decay spectra is of the order 



MJ 



V 



< 1 



(4.6) 



where mi is the mass of the final charged lepton. The other form factors are 



2P+/1 
2P+K^f2 

2P^gi 
2P+K^g2 



Mi(B',]\v+ B,i 
-M,(5',T A+ B,i 



(4.7) 



What is usually measured is the total decay rate R, the electron-neutrino correlation 
aeu and the electron Og, neutrino ai, and final baryon as asymmetries. The e-u correlation 
is defined as 



"^'^-'iV(0,,<i7r) + iV(ee.>|vr) ' ^^-'^ 

where N[Qeu < l^r) is the number of e-u pairs that form an angle Qeu smaller than 90°. 
The correlations ae,au and are defined analogously with Qe,Qu and 9^ now being the 
angles between the e, z/, B directions and the polarization of the initial baryon. 

Ignoring the lepton-mass one can calculate expressions for the measured quantities. 
We copy the expressions for R, aev, cte, Oiu and as from Ref. |]57|: 



ei/ < 2^) 



iv(e 



R 



G' 



AM5|\/|2 
607r3 



lp + \p')fl + ^^P'f2+{^-lp + yP')9l 



Rap 



Rap 



+Y(3 92 + 7/571/2 + (-4/? + Qngig2 + (/lA/ + 5giXg) 



G 



607r3 



;i-^/5 + y/5^)/i^-^/5V2^ + (-l 



2 24 
-2P'9l - ^P'fif2 + (4/? - 2p')g,g2 - -yP'diK 



(4.9) 



(4.10) 



G 



4 



Ip + ^^P')fi - Yi^'fl + (-2 + ^^P')9l - -/9I 



607r3 

+ (-^/5 + Yl^')f29i + ^/5V2^72 + (y/5 - |/5^)^7l^72 
4 



+-/?'(^iA/ + hXg - Ag^\ 



(4.11) 
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Ray 



a 



2AM5|\/|2 



607r3 



+ (3/5 - )/l/2 + (2 - + -nhgi + (-3/5 + )/l^72 

+^/3'((7iA/ + /iA3 + 4(7iA, 



Ra 



B 



G- 



6 



3 o 3 3 



(4.12) 



(4.13) 



where p is defined as /3 = (M^ - M/) /M^, and AM = Mi- Mf, Mi, Mj being the masses 
of the initial and final baryon, respectively. The dependence of /2 and 92 is ignored 
and /i and (71 are expanded as 



(4.14) 



We get the corresponding expression for the dipole parameterization by putting 

\f = 2M^f,/M^, \g = 2Mfg,/Ml. (4.15) 



Two corrections have to be made to these quantities: (a) the correction due to the 
nonvanishing lepton mass and (b) the radiative corrections. To keep the lepton mass 
non-zero we integrate the differential decay rate for f2=92 = ^'- 



R = G 



,Ml\Vl 
327r3 



{l-x'f 



dz 



{fl + 9t)h+{fl-9l)h 



(4.16) 



with 



h 
h 
f 

X{x,y,z) 



— {z — Xi)x f , 

X'/'il,z,x")\'/\z,xlO)/z, 
x^ + y'^ + z^ — 2xy — 2yz — 2xz 



mi 
Mi' 



X 



Ml 
Mi 



= zM? 



The ratio R/R{mi = 0) for the various reactions is given in Table |4.1| for the e-mode 
(/ = e) and the /i-mode (/ = /i). 

The radiative corrections are well known |57|. The rate can be written as 

R! = R{l + 5a){l + 5b) = R{l + 5) , (4.17) 

where R is the rate defined in Eq. (|4.9|) . The term 5a comes from the model independent 
virtual corrections and the bremsstrahlung pS. For the model dependent term 5;, one 
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Table 4.1: Ratio of the rate to the rate with vanishing lepton mass. 



rieaction 


e- 


R/R{mi = 0) 
-mode /x-mode 


np 





.472 




E+A 


1 


.000 




S-A 


1 


.000 




E-EO 





.955 




EOE+ 





.830 









.971 




Ap 


1 


.000 


0.161 


eV 


1, 


.000 


0.431 


E^n 


1 


.000 


0.443 


S-A 


1 


.000 


0.271 


S-EO 


1 


.000 


0.0136 


50E+ 


1 


.000 


0.00821 



Table 4.2: Radiative corrections to the semileptonic decay rates. 



Reaction 




6 


e-mode np 


0.0486 


0.0706 


E+A 


0.0015 


0.0225 


E-A 


0.0012 


0.0222 


E"E+ 


0.0226 


0.0441 




0.0104 


0.0316 


Ap 


0.0207 


0.0421 


TPp 


0.0196 


0.0410 


Ti~n 


-0.0025 


0.0184 


S-A 


-0.0015 


0.0195 


S-EO 


0.0000 


0.0210 


/x-mode Ap 


0.0468 


0.0688 


TPp 


0.0336 


0.0553 


E-n 


-0.0022 


0.0188 


S-A 


-0.0013 


0.0197 


S-E° 


0.0002 


0.0212 
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takes the value 0.021 [^. The 6 term is the whole radiative correction given in Table ^ 
together with 6a. 

For the angular correlation aeu and the asymmetries ae and ai, the model dependent 
corrections vanish and the model independent corrections are of the order of 0.001 |5T 
We therefore do not include these corrections. 
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Table 4.3: Matrix elements for weak beta decay. 



Reaction 

-({x''\03\x''} + {x''\03\x''}) 
S+A Ux'' IO3I x'') + (x'' IO3I x'')) 

S-A f (x^i IO3I x'') + [x'^ m X'V) 

Ap ^Ax^^-x^^mx'"") 

75(x'' + x''|03|x'') 

2-A ^((x^^|03|x'V(^''l^3|x'')) 

-;^((x^M03lx^^) + (x^^|03|x^i)) 

2°s+ -{{x''mx'') + {x''mx'')) 



4.2 The form factors in the quark model 

The basic formula for the matrix elements in Eq. (^4.7|) is Eq. ( p.25| ). The Dirac quark 
current 

qV^q, V^ = ril-1^) (4.18) 
can be generalized by the Ansatz 

r'^ = h,l^ - ^la'^'K, + ^K'^ + gi.ri, - ^-^la^'K,^, + ^-^K'^^, . (4.19) 
niq nig nig nig 

The subscript 'g' stands for a transition on the quark level. The form factor f2q is de- 
termined by the anomalous quark moments through CVC But f2q as well as f^q, g2q 
and 5^35 do not contribute to /i(0) and 5'i(0) because of their factor K. Since K'^ is 
small and contributions different from /i(0) and gi{0) are of higher order in j3 we put 
f2q = f3q = Qiq = dsq = wlthout loss of generality. For the form factors fig and gig we 
shall distinguish between the transitions d ^ u and s ^ u. In the last Chapter we have 
put the parameters fiud and giud equal to one. For the s ^ u transition we note that we 
have to put fius = 1, if we wish to predict \4s, since the rate is proportional to /i^^lVusp. 
The determination of gius is discussed in the next Section. Notice that we would have to 
put fius ~ gius ~ 5, if we liked to fit the weak sector with parameter set 4 (see Table 



As seen from Eq. ( [4.7|) we have to calculate (x| jV^"^] Xt), (xt I^^IXi)? (Xt I^'''Ixt) 
and (xtl^^lXi)- Table O gives the matrix elements for the various reactions with 
0,, = V+,A+. 
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Table 4.4: Parameters in Eq. ( [4.201 ). 



Reaction 


Aifi) 


Aif2) 


Ai9i) 


np 


1 


{2A2 - 5Ai)/3 


5 

3 


S+A 





{A2 + A3-2A,)/VQ 


[2 

V 3 


S-A 





{A2 + A;-2A,)/VQ 


[2 
V 3 




V2 


-{AA, + A2 + As)/{3V2) 


3 


EOE+ 


-V2 


iAA, + A2 + As)/i3V2) 


2v^ 




-1 


(2A2 + 2A3- Ai)/3 


1-^ 
3 



For = we have for AS* = transitions 
/i = ^(/i), 



/2 



(4.20) 



(2 



62 + Qi ' 



with As given in Table [4.4| . The values A{fi) and ^(s'l) are the values in the nonrelativistic 
quark model. 

The AS" = 1 transitions for = are 



/i 



{2nY 



d^d-'Q 



EW±\ $t(MO$(M)i?(/i) 



(4.21) 



{2ny 



E'sE[,M \ $t(M0$(M)i?((7i) 



5(/i) = 5i(a'a + Qi)2(6'6 + gi) 

+i?2(a'-a)2Ql(6'6 + Q2^ M - 



c2 + gi)(rf2 + g2) 

/ 2 ^2 

+53(a' - a)(6' - b)Ql{a'a + Ql) 



2 / ' 



(4.22) 



B{gi) = B,{b'b~Q'_ 



ia'a + Ql)^ + ia' -afQ^ 



{cd-qlf 



(c2 + gi)(d2 + gi) 
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Table 4.5: Parameters in Eq. ( [4.221 ). 



Reaction Bi B2 -B3 B4 Bq 

vo„ _J_ _±_ V2 _±_ WI V2 

^ sy^ 3 3^2 3 3 



S-n -1 



_ 2 1 8 2 

3 3 3 3 3 

'=:- \ /3 J_ L J_ _9 /2 _i 

" V2 v^v^ y 3 6 

"-VO J_ _5_ _J_ _5_ 

V2 3V2 3\/2 3\/2 3\/2 6 

5 15 4 1 

_3 3 3 3 3_ 



1 



+55(a'-a)^Qi(6'6-Q 



cdqj 



+BQ{a' -a)Ql{b' + b){a'a + Q 



c2 + gi)(rf2 + g2) 

/ -2 ^2 



+ 



The Bi for the different decays are given in Table 



Eqs. ( [4.21| ) and (|4.22|) confirm the Ademollo-Gatto theorem |6^. Since [a' — a) ~ Am 
and {b' — 6) ~ Am the symmetry breaking for fi is of the order (Am)^ whereas it is of 
the order Am for gi owing to the term containing Bq. In addition to Ademollo-Gatto we 
see that the symmetry breaking for gi{A p) is also of second order. 

The full formulae for < are longer than the ones for = 0; they are given in 
Appendix ??. 
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4.3 Results 



The form factors can be determined by the generahzation of Eqs. ([4.20|) and ( [4.21 
the parameterization of the form factor f{K'^): 



m 



l-XyAj + K^/A 



With 



(4.23) 



we get the result shown in Tables [4.6| and together with the rates, angular correlation 
and asymmetries from Eqs. ( [4.9|) -( ^rT^ ). The parameters A„ are determined by the calcu- 
lation of the appropriate derivatives of fiK"^) at = 0. The rates have been corrected 
taking into account the non-vanishing lepton mass and radiative corrections. 



4.3.1 The rates and /i(0),^i(0) 



We fit our remaining parameters rris and /3s = /3a to the rate R and the ratio gi/fi for 
the processes A — > pe~Ue and —>■ ne~Ue, and /3h is fitted to the rate for Ae~Ue- 



We find parameter set 5 of Table 3.1 and get 



R{A 
i?(S~ 

9i/m 

9i/fi{^- 



pe Ve) 

ne~iJe) 



3.38 X lO^'s" 
5.79 X lO^s" 
0.782 , 
-0.261 . 



But the value for seems to be too small if we compare it with the well confined 
value in the meson sector of the same model . By considering the constraints Am = 
— ^u/d — 140 MeV and ms/rriu/d — 1-4 @ we choose = 0.40 GeV (set 6). The 
results with set 6 have been collected in Tables Bl6| and ffl7. 



The largest discrepancy between theory and experiments comes from the rates and 
Qil fi for the processes A — >■ pe~Vf, and S~ ne~Ue- By changing the axial couplings of 
the quarks, i.e. gius — 0.9, we could improve the rates of both reactions, but the ratios 
(71 //i clearly force us to use gius = 1- Another modification could be the A — E°-mixing. 
Let us write 



A 



phys 



AA + BE\ Sjhy. 



-5 A + A S'- 



A^ + B^ = 1. 



(4.24) 



From the measurement /i/(y'i(E — * Ae Ue) we get a constraint A > 0.9961, B < 0.0078 
which gives 



fi/gi{A pe-Ue) > 0.773 (> 0.73 for Am = 63 MeV) 



(4.25) 



and reduces the rate by about 1 %. Therefore, A — S°-mixing only improves one of the 
four values for which theory and experiment differ. Actually, this inconsistency of our 
values is a general feature of quark models with a SU(6) flavor-spin symmetry.^] The ratio 



^The bag model calculation 
ne've) = -0.28. 



gives similar results: gi//i(A 



0.84, and 
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Table 4.6: Results for AS = weak beta decay with parameter set 6. Experimental data 



are from PDG 







np 


S+A 


S-A 








fl 


/i(0) 
Ai [GeV] 
As [GeV] 


1.00 
0.69 
0.96 




-0.32^ 
-1.72^ 




-0.32^ 
-1.72^ 


1.41 
0.60 
0.81 


-1.41 
0.60 
0.81 


-1.00 
0.56 
0.71 


91 


51 (0) 
Ai [GeV] 
Aa [GeV] 


1.25 
0.76 
1.04 


0.60 
0.77 
1.05 


0.60 
0.77 
1.05 


0.69 
0.77 
1.04 


-0.69 
0.77 
1.04 


0.24 
0.76 
1.04 


Qi/h 


Theor. 
Expt. 


1.252 
1.261 
±0.004 


0.736'^ 
0.7421^ 
±0.018 


0.736^ 


0.491 


0.491 


-0.244 
< 2 X 10^ 




Theor. 
CVC 


1.81 
1.85 


1.04 
1.17 


1.04 
1.17 


0.76 
0.60 


-0.76 
-0.60 


0.73 
1.00 


f [GeV-i] 






















Rate [lO^s-i] 
e-mode 


Theor. 
Expt. 


1.152 X 10-9 
1.125x10-9 
±0.004 


0.24 
0.25 
±0.06 


0.389 
0.387 
±0.018 


1.47 " 


3.65^ 


1.55^= 





Theor. 


-0.101 


-0.404 


-0.412 




Expt. 


-0.102 


-0.35 


-0.404 






±0.005 


±0.15 


±0.044 




Theor. 


-0.112 


-0.701 


-0.704 




Expt. 


-0.083 










±0.002 








Theor. 


0.989 


0.647 


0.645 




Expt. 


0.998 










±0.025 






as 


Theor. 


-0.548 


0.070 


0.077 



0.438 



0.288 



0.793 



-0.514 



Expt. 



Instead of A, we list ^ Instead of gi/ fi we Ust y^^i. " xlO'^ ^ xlO"^. 
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Table 4.7: Results for AS = 1 weak beta decay with parameter set 6. Experimental data 



are from PDG H 







Ap 




S n 








fl 


/i(0) 
Ai [GeV] 
A2[GeV] 


-1.19 
0.71 
0.98 


-0.69 
0.64 
0.84 


-0.97 
0.64 
0.90 


1.19 
0.68 
0.89 


0.69 
0.75 
1.05 


0.98 
0.75 
1.05 


91 


9m 

Ai [GeV] 
As [GeV] 


-0.99 
0.81 
1.12 


0.19 
0.83 
1.16 


0.27 
0.83 
1.16 


0.33 
0.81 
1.10 


0.94 
0.81 
1.12 


1.33 
0.81 
1.12 


9\lf\ 


Theor. 
Expt. 


0.826 
0.718 
±0.015 


-0.275 


-0.275 
-0.340 
±0.017 


0.272 
0.250 
±0.042 


1.362 
1.287 
±0.158 


1.362 
< 2.93 




Theor. 
CVC 


-0.85 
-1.19 


0.44 


0.62 
1.12 


0.070 
-0.080 


0.98 
1.38 


1.38 
1.95 


§ [GeV-i] 




-0.025 


0.0043 


0.0061 




_a 


_a 


Rate [lO^s-i] 
e-mode 


Theor. 
Expt. 


3.51 
3.170 
±0.058 


2.72 


5.74 
6.88 
±0.26 


2.96 
3.36 
±0.19 


0.549 
0.53 
±0.10 


0.942 


Rate [lOS-i] 
/Lt-mode 


Theor. 
Expt. 


0.58 
0.60 
±0.13 


1.18 


2.54 
3.04 
±0.27 


0.80 
2.1 
±2.1 


7.47 X 10^3 


7.74 X 10~3 





Theor. 


-0.100 


0.443 


0.437 


0.531 


-0.252 


-0.248 




Expt. 


-0.017"^ 




0.279 


0.53 










±0.023 




±0.026 


±0.1 








Theor. 


-0.021 


-0.536 


-0.537 


0.236 


-0.226 


-0.223 




Expt. 


0.125 




-0.519^= 












±0.066 




±0.104 








ay 


Theor. 


0.992 


-0.318 


-0.318 


0.592 


0.973 


0.973 




Expt. 


0.821 




-0.230^ 












±0.066 




±0.061 










Theor. 


-0.582 


0.568 


0.569 


-0.519 


-0.437 


-0.439 




Expt. 


-0.508 




0.509= 












±0.065 




±0.102 









a 92 
91 M 



0.023 since ^ 

91 



constant. ^ From Ref. IGTI]. From Ref. 
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Qi/fi can generally be written as 



91 
fi 




(4.26) 



non— rel 



where (fi'i//i)non-rei is the non-relativistic value. The quantity p is a relativistic suppres- 
sion factor due to the " small " components in the quark spinors (in the bag-model) or 
due to the Melosh-transformation (in our model). The quantity r] is an enhancing factor 
due to SU(3) symmetry breaking in AS = 1 transitions. From Tables ^]6| and |4.71 we 
see that p ~ 0.73 — 0.76 depending on the strangeness content of the wave functions and 
?7 ^ 1.11. This simple estimate shows that every quark model is a priori constrained to 



^i//i(A->pe Ue) ^ _g 
5'i//i(S" ne-Pe) 
in contrast to the experimental value —2.11 ±0.15 for g2 



that 



91 
fi 



0.715 + 0.28 



Ap 



92 

fl 



(4.27) 

0. For 7^ it is measured 

(4.28) 



and 62 



91 _ Q_237.^2 



/: 

which will bring the data closer to 
small to remove the discrepancy. 



fl 



0.34 ±0.017 



(4.29) 



-3, but in our model 5'2/fi'i — 0.025 which is much too 



4.3.2 /2(0) and ^2(0) 

Our model agrees with the conserved vector current (CVC) hypothesis. The deviations 
have the same origin as the too small neutron magnetic moment since /2 and the magnetic 
moments have similar analytic forms. If we take fip as in Sec. ^.2| , CVC will reproduce 
our values. The experimental situation is not yet clear, some experiments confirm 
and some disprove ||6T| CVC. 

For AS" = 1 transitions the prediction of g2/9i for nonrelativistic quark models is 
~ 0.37 and for the bag model ~ 0.15 [^. Our model gives also a constant value 



92 

.91, 



AS=1 



For AS" = transitions we get 




0.025 



0.0033 



(4.30) 



(4.31) 



if we put md — mu = 7 MeV. This confirms the viewpoint of the PDG ||5^ which fixes 
92 = 0- Experiments also find a vanishing or small g2 |57[ . 

With CVC and the absence of g2 we reach the same conclusion that was reached in 
nuclear physics. 

^(ff2)HA and (32)51; could only be calculated for Am = 63 MeV and we get {92/91)^3=1 — 0.11. 
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1.2 



I 0.8 

t 

C3 



M 0.4 




Figure 4.1: The axial vector form factor gi{K'^) for the np-transition. The dipole formula 



is compared with the experimental data taken from Ref. [35 



4.3.3 ii'^-dependence of the form factors 



Tables ^]6] and suggest that the form factor of Eq. ( [4.23| ) can be approximated by the 
dipole form 



/(O) 



(4.32) 



(1 - KyAif ■ 

The axial vector form factor gi for the neutron decay gives a value Ma = A2 = 1.04 
GeV compared to the experimental value = (1.00 ± 0.04) GeV Figure [4.1| 

compares the dipole-fit with experimental points ||65|| . 

If we take the dipole Ansatz we can compare our values for My and with the 
results of other work (see Table |4.8| ). 

The contribution of My and Ma to the rate and to x = gi/ fi to first order is 



AR 

~r" 

Ax^ 

x"^ 



8 



7(1 + 3x2) \M, 



l_ ^ 

^ Ml 



V 



(4.33) 



7 



1 - Ueu) 6 + 5 



M2 



Ml 



which shows that our parameters give for the decay — ^ ne^Ue a 0.3% larger rate and 
a 5% smaller gi/ fi than with the parameters of Gaillard et al. that are often used for 
the experimental analysis. Although this does not explain the inconsistency of the data 
with our calculation, it shows that future high-statistics experiments should pay more 
attention to My and Ma in analyzing gi/fi- 
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Table 4.8: The parameters My and Ma for varfous models. 



This work Gaillard et al.|59| Garcia et al. |57| Gensini [54 





My 


Ma 


My 


Ma 


My 


Ma 


My 


Ma 


np 


0.96 


1.04 


0.84 


1.08 


0.84 


0.96 


0.84 


1.08 


SA 




1.05 




1.08 




0.96 




1.08 


EE 


0.81 


1.04 


0.84 


1.08 


0.84 


0.96 


0.84 


1.08 




0.71 


1.04 


0.84 


1.08 


0.84 


0.96 


0.84 


1.08 


Ap 


0.98 


1.12 


0.98 


1.25 


0.97 


1.11 


0.94 


1.16 




0.84 


1.16 


0.98 


1.25 


0.97 


1.11 


0.94 


1.16 


En 


0.90 


1.16 


0.98 


1.25 


0.97 


1.11 


0.94 


1.16 


HA 


0.89 


1.10 


0.98 


1.25 


0.97 


1.11 


0.94 


1.16 


HE 


1.05 


1.12 


0.98 


1.25 


0.97 


1.11 


0.94 


1.16 



4.4 Cabibbo fit and K 



us 



There are at present some questions concerning flavor SU(3) breaking in semileptonic 
weak hyperon decay 0, 0, 54 1. This symmetry breaking is included to all orders in our 



approach. We can extract this effect from our model and fit the experimental values 
within the Cabibbo model. This has been done for the bag model in 1987 but 
since then some of the experimental data have changed. The main difficulty is the rate 
for S~ —* Ae~Ue, which comes out too small and is made even worse by their symmetry 
breaking scheme. 



The reason for using the Cabibbo model is the bad fit for the relevant data in Table ^7? . 
The experimental deviations are as large as 17% and do not permit a precise determination 
of the Kobayashi-Maskawa matrix element Vus- 

The Cabibbo model analyses the experimental results in terms of three parameters 
Vus, F, and D. The form factor /i(0) is given by SU(3) symmetry, /2(0) by its CVC-values 
and gi{0) as follows: 

np Ap Tji S^A ha H^S 

F + D -,J3/2{F + D/3) -F + D ^J2/3D ,J3/2{F - D/3) 1/V2{F + D) 

The purpose of a Cabibbo fit is to determine the values of Vus, F, and D corresponding 
to the best agreement between experiment and theory. 

We choose the value Vud = 0.9735 of Ref. |Q, and perform a least-squares x^-fit. 
Radiative corrections are taken into account. 

Table ^]9| gives the residuals Ri = [xj(fit) — Xj(meas)] /(5xj(meas) for each reaction, 

= J2i Ri for 11 degrees of freedom, and the fit variables. 

Two fits have been made without symmetry breaking. Fit 1 uses the commonly used 
values of My/A given in Ref. and fit 2 uses our masses, which gives a slightly improved 



X^. Therefore, we use our values for all the other fits. The rate for S —>■ Ae u produces 
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Table 4.9: Fits to the Cabibbo modeL The decay rates are given in units of 
10^s~^, except for the neutron decay with 10^s~^. The residuals are given by 
[xj(fit) — Xj(meas)] /5xj(meas). 



Expt. 




Residuals of fit no. 








1 


2 3 4 5 


6 


7 



R{n - 
i?(A - 
i?(S - 

RiE- 
R{E- 
R{A - 
i?(S - 
RiE- 
9i{n - 



pep) 1.125 ±0.013 -1.1 -1.1 -1.2 -1.5 -1.2 -1.2 -1.1 

pep) 3.170 ±0.058 0.9 0.8 1.0 1.5 1.2 0.3 0.4 

■nev) 6.88 ±0.26 -0.5 -0.3 -0.3 -1.1 -0.5 -0.5 -0.7 

-^Aei/) 0.387 ±0.018 2.7 2.7 1.1 2.2 0.4 1.2 1.1 

Aeu) 0.25 ±0.06 0.2 0.2 -0.1 0.1 -0.2 -0.1 -0.1 

-^Aez/) 3.36 ±0.19 -2.2 -2.1 -2.8 -3.5 -3.3 l.P 1.1^ 

-^Sez/) 0.53 ±0.10 0.0 -0.1 -0.2 0.2 -0.3 -0.3 -0.3 

p/ii/) 0.60 ±0.13 -0.3 -0.3 -0.3 -0.3 -0.3 -0.4 -0.4 

•n/ii/) 3.04 ±0.27 -0.2 -0.1 -0.1 -0.4 -0.2 -0.2 -0.3 

-^A/iz/) 2.1 ±2.1 -0.6 -0.6 -0.6 -0.6 -0.6 -0.6 -0.6 

■pep) 1.261 ±0.004 -0.1 -0.1 -0.2 -0.9 -0.3 -0.1 0.1 

>pep) 0.718 ±0.015 0.9 1.0 2.8 5.2 3.5 2.7 1.1^ 

>nev) -0.340 ±0.017 0.9 1.0 -0.6 -1.8 -1.4 -0.8 -0.5 

Aeu) 0.25 ±0.05 -0.9 -0.9 -1.1 -0.8 -1.1 -1.1 -1.1 



F 
D 

Ks [±0.003] 



0.468 
0.793 
0.227 
17.5 



0.469 
0.792 
0.226 
16.9 



0.461 
0.799 
0.226 
21.4 



0.458 
0.800 
0.222 
54.8 



0.457 
0.802 
0.223 
30.3 



0.460 
0.801 
0.225 
14.3 



0.462 
0.799 
0.225 
7.7 



^We take the experimental value 1.83 ± 0.79 
taken. 



69 . ''The effective value 0.744 ± 0.015 is 



Table 4.10: Symmetry breaking for /i. The ratio fi/jf^^^^ is shown 



Am 



This work 
63MeV Am -- 



Donoghue Krause 



133MeV 



AS = 
Ap 

Era 
EA 



1 

0.978 
0.979 
0.978 
0.981 
0.982 



1 

0.976 
0.975 
0.975 
0.976 
0.976 



1 

0.987 
0.987 
0.987 
0.987 
0.987 



1 

0.943 

0.987 
0.957 
0.943 
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Table 4.11: Symmetry breaking for gi. The ratio Qi/g^^^^ is shown 





This work 
Am = 6oMeV Am 


= IddiVleV 


Donoghue 


np 


1.000 


1.000 


1.000 


SA 


0.959 


0.981 


0.9383/0.9390 


SS 


0.955 


0.982 






0.916 


0.977 




Ap 


1.021 


1.072 


1.050 




1.011 


1.051 




Sn 


1.012 


1.056 


1.040 


SA 


0.987 


1.072 


1.003 


SE 


0.981 


1.061 


0.9954 



the largest deviation, indicating SU(3) symmetry breaking. Using Tables 4.10| and 4.11 



we fit using Am = 63 MeV (fit 3), Am = 133 MeV (fit 4), and the results of Donoghue 
(fit 5). The symmetry breaking scheme for Am = 133 MeV cannot be correct. We 
improved the above deviation at the cost of introducing a large contribution to from 
the rate H~ —>■ Ae^u and the ratio gi/fi for A —>■ pe'v. Fit 3 performs best, but is 
still worse than no symmetry breaking. Since there are doubts about the experimental 
rate for S" Ae'v we use (1.83 ± 0.79) x lOS"^ H for fit 6 with Am = 63 MeV. 



This case is nearly as good as no symmetry breaking, which gives = 13.8. We have to 
remember Eq. ([4.28|) and footnote on page and the nonvanishing g2 gives an effective 



gi/fi = 0.744 ± 0.015. With these data, fit 7 is in excellent agreement with experiment 
(x^ = 7.7 for 11 DF). We observe that Vus lowers its value from fit 1 to 7, so that we get 
closer to the one derived from the meson sector. The analysis of i^es decays yields 

^ r 0.2196 ±0.0023 (Ref. [32]) , 

1 0.2199 ±0.0017 (Refs. 0). ^ ^ 

Other values from hyperon beta decay are 

_r 0.2258 ±0.0027 (Ref. [4]) , 

1 0.222 ± 0.003 (Ref. 0) . ^ ' 



The last value is derived from WA2 data |71| with the symmetry breaking scheme from 



Ref. [^|. Our result is in excellent agreement with the recent value from Ref. and it 



is also consistent with the other baryonic value [52|. But the discrepancy with the meson 
sector still remains. 

In conclusion, there is evidence for symmetry breaking from the rate E~ Ae~z/. 
But symmetry breaking alone makes the fit even worse. We also have to use a new value 
of My/A and a small second-class axial coupling as given by our model. In this Section, 
we considered symmetry breaking due to the mass difference m^ — rriu/d alone, but an 
asymmetric wave function as treated in Chapter ^ breaks SU(3) symmetry as well. 
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4.5 Conclusions 

Parameter set 6 yields the best fit for the semileptonic weak decays in the MQM. It agrees 
with the data within 17%, except for some of the correlations and asymmetries. The main 
discrepancy with experiments lies with the decays A — > pe~i'e, and E~ ne~i>e, which 
are the main sources for determining the Kobayashi-Maskawa mixing matrix element T4s- 
This is the reason why we switched over to the Cabibbo model in calculating Vus- 

The parameter set 6 is quite different from the set 4 found by a comparison with 
the baryon magnetic moments. This result shows that the representation of the model 
discussed up to now does not give a consistent picture of the u, d, and s sector of the 
baryons. A consistent picture can be achieved with the help of asymmetric wave functions 
as shown in the next Chapter. Even the fit within the weak sector improves dramatically. 
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Chapter 5 

Asymmetric wave functions 



5.1 The diquark model 

In this Chapter we investigate the effects of an asymmetric wave function. Since the early 
days of the quark model diquark clustering has been studied. In Gell-Mann's original 



paper on quarks he mentions the term diquark in a footnote. Refs. |72, 73 1 first took 



the idea seriously, and in both papers, a baryon is described as a bound state of a quark 



and a diquark. A recent review with many references can be found in Ref. [Q. It has 
even been seen that many gluon effects can be simulated by diquarks [|75| . 



The concept of diquarks is also useful in treating deep-inelastic electron scattering. 



Feynman []T6[ observed that the experimental ratio of the neutron to proton structure 
function can be qualitatively explained if nearly all the momentum is carried by a leading 
quark, which is a m quark in the case of a proton and a d quark in the case of a neutron. In 
both scalar diquark with isospin zero remains. Close includes both scalar and 

vector diquarks, which have the same probability within SU(6) symmetry. The conclusion 
from the analysis |]7^ is that a scalar diquark is more probable than a vector one at large 
momentum transfer. This is in agreement with the fact that the QCD spin-spin force, 
originally introduced into the quark model in Ref. |]TS[, is attractive and strongest in the 
spin-0 quark-quark state. 

Most generally the proton can be composed of the four combinations for spin and 
isospin either equal to zero or one: 

\p) = A u{ud)Q(f) + B d{uu)o(j) + C u{ud)i(f) + D d{uu)i(f) + perm., (5.1) 
+ + C^ + = 1 , 

where the parentheses indicate the diquark clustering and the spin of the quark pair is 
given as a subscript. A preliminary calculation for the magnetic moments of the baryon 
octet suggests that C and D are small or zero. The fit shows that if we put B = -0.2 the 
parameters Pg and Pq, discussed below, can be chosen to be equal. Note that SU(6) is 
still broken in this case. 

Considering the number of the parameters and because of the picture from deep- 
inelastic electron scattering we only keep the spin-isospin-0 part. To implement this 



feature we have to modify the wave functions in Section 2.3. In order to get a spin 
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isospin-0 clustering we write the proton wave function as 

-1 



iX'' - hx'') + duu + hx'')] ■ (5.2) 



The wave functions for the Ss and Ss are obtained by changing the flavor wave function 
accordingly. The A wave function is given by 



-1 



(paX'^^ i'^ds — dus) + 02X^^ {usd — dsu) + (pix''^ {sud — sdu) 



(5.3) 



The function 0j is the momentum wave function with symmetry in the quarks different 
from i. We choose 

(t>i = 7V,e-^* , (5.4) 

where the normalization factor Ni is given below and the Xi are the generalized forms of 
M72/3: 



^3 
X2 



+ 



m? 



+ 



2 ' 



2/32/32^^(1- 0(1 -r^ + e^) ^^2i5li5lr){\-ri){\-r)-ri^) 



+q±Q. 



PI -PI 



2 

mi 



ml 



+ 



mt 



Pl^P^^r]{l - ri + Crj) 2r^i(5l 277(1 - 0/3§ 2(1 - 77)/?; 



^1 



9± 



2 (1 - OPl + ^(1 - riWq + g2 (1 - 0/?^ + ^(1 - 



2/32/32^^(1- 0(1 -^r7) 2/3^/3277(1 -77)(l-er/) 



32 _ 02 



iir- 



111: 



Wvi'^-^V) 2r7^/^Q 277(1 -0/9,^ 2(1 - 77)/32 
There is a special form of iVj with Ni = N2 = N^: 



(5.5) 



iV/ 



2X3 d$,dq±drjdQ± 
^77(1- 0(1-^) 



(5.6) 



However, this normalization gives too small rates for the semileptonic decays. We use 
instead our usual normalization 



j d\d^Q\cPi 



1 . 



(5.7) 



5.2. RESULTS AND DISCUSSIONS 



51 



5.2 Results and discussions 

The calculations of the matrix elements are similar to those in the symmetric case, but 
they are more involved because the wave function is no longer symmetric. We give an 
example for the vector current [K"^ = 0) [see Eq. ??] 



(tn I TTi) asym (tn I TTi) sym 

2(a - a'fibb' + Q±)cd{qLQBf 



a'2 + Ql){a^ + Ql)Jb'^ + QlJh^ + Ql{c^ + ql){d? + ql) 



(5i 



For the numerical calculation it is important to simplify the additional term to reduce 
the six-dimensional integral to a five- dimensional one (see Appendix [A.2| ). 



Before presenting the results we give some general considerations. To take an asym- 
metric wave function is the most natural extension to the minimal model. Since we do not 
have this possibility in the meson sector, this extension does not contradict the minimal 
meson model ||10|, [1^. Another question concerns the ratio in Eq. (|4.27|) . It is in principle 



possible to improve this value drastically to give —2.26 {rriu/d = 0.26 GeV, = 0.39 
GeV, I3qn = 0.3 GeV, /S^n = 0.7 GeV, /3qs = 0.7 GeV, /3g^ = 0.3 GeV), which shows 
great flexibility compared to other models. 

We minimize the function A = I]j |xj(fit) — a;j(meas)| for the fit. The details of 
the fit are described in Appendix |A.3| . Figure ^.1| shows A as a function of the various 
parameters. The black areas in the density plots indicate the minimum, the white areas 
the maximum of the function A. The fixed values of the parameters in the plots are the 



ones in Table f).l\ The figures (a) and (b) show the fit for the magnetic moments of the 
proton and neutron, and the weak axial form factor gi{n pe'v^)- We can see that 
there is no large diquark clustering in the nucleon sector (/Jqat = l^qN)- The magnetic 
moments of the hyperons are fitted in figures (c), (d) and (e). There is a strong diquark 
clustering in the strange baryon sector {jSq ~ 2/3q). The mass of the u and d quarks is 
(0.26 ± 0.20) GeV for both the nucleons and hyperons. The strange quark mass is fixed 
to be 1.5 times niu/d- Figure (f) gives the fit to the rates and the axial form factors for 
the decays A pe~Ue and E ne~Ue as a function of the weak quark form factors. 
We summarize the results for the asymmetric parameter set in Table There is a 



considerable improvement for the magnetic moments of the neutron, the Ss and Hs, and 
for the rates and the ratios gi/fi of the semileptonic decays. Non-zero quark anomalous 
magnetic form factors could even give better results for the magnetic moments of the 



neutron and the H (see Table 3.3). Parameter set 8 fits both electromagnetic and weak 



Table 5.1: The parameters of the asymmetric constituent quark model. All numbers are 
given in GeV. 



"fT^u/d f^s PqN PqN PqS/A /^gS/A PqB PqE. flus dlus 



Set 8 0.26 0.395 0.55 0.55 0.35 0.65 0.34 0.65 1.32 1.17 
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Figure 5.1: Parameter space of the asymmetric wave function fit. The deviation from 
the experimental data is plotted against the various parameters. The black areas in the 
density plot show the minimum, the white areas the maximum of the difference between 
the experimental values and those given by the fit. The /3s and the masses are given in 
units of GeV.(a), (b) nucleon sector; (c), (d), (e) magnetic moments of the hyperons; (f) 
semileptonic, strangeness changing, weak decay. 
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Table 5.2: Electroweak properties of the baryon octet. The calculations with parameter 
sets 4, 6 and 8 are compared. We see that set 4 is only able to fit the magnetic moments 
(see Section p?2|) , set 6 fits the weak sector (see Chapter and set 8 can fit both sectors 
simultaneously. The magnetic moments are given in units of the nuclear magneton, the 
decay rates are given in 10^s~^. 



Particle Experiment |^| Set 4 Set 6 Set 8 



/.(E+) 

/i(A) 

/.(S-) 

t(Ap) 
l(Srz) 

t(H-A) 

^(S-SO) 

R{Ap) 

i?(En) 

R{E-A) 

i?(S-SO) 



2.79 ±10-^ 2.85 2.78 2.82 

-1.91 ±10"^ -1.83 -1.62 -1.66 

2.42 ± 0.05 2.59 3.23 2.48 

-1.157 ± 0.025 -1.30 -1.36 -1.09 

-0.613 ± 0.004 -0.48 -0.72 -0.64 

-1.250 ± 0.014 -1.25 -1.87 -1.28 

-0.679 ± 0.031 -0.99 -0.96 -0.78 

1.261 ± 0.004 1.63 1.252 1.25 

0.718 ± 0.015 0.957 0.826 0.760 

-0.340 ± 0.017 -0.319 -0.275 -0.238 

0.250 ± 0.042 0.319 0.272 0.190 

1.287 ± 0.158 1.594 1.362 1.13 

3.17 ± 0.058 0.14 3.51 3.20 

6.88 ± 0.26 0.16 5.74 6.68 

3.36 ± 0.19 0.10 2.96 3.61 

0.53 ± 0.10 0.02 0.55 0.48 



properties of the baryon octet. Weak form factors of the quarks are needed for the 
strangeness changing transition s — >• m. Since the rate is proportional to \ Vus\^fius we 
have to fit both values, we cannot predict an accurate value for Vug. 
General features exhibited by these results are 

• In the nucleon sector, there is no diquark clustering {Pqn = Pqn)- 

• There is a strong diquark clustering in the strange baryon sector (2/?q ~ Pg). 

• The momentum scale parameter for the diquark pair is about the same for all 
baryons (Pqn = PgN ~ PgS/A = Pqe)- 

Table ^.3| compares the results of this work with other models. The nonrelativistic 
quark model yields accurate magnetic moments, but fails in the weak decay sector. The 
results derived from QCD sum rules and the bag model are comparable to those from our 
model. Nevertheless, some of the data are reproduced best within our model. The fits 
from the lattice simulation and the Skyrme model are too small. 
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Table 5.3: Electroweak properties of the baryon octet. The calculations of the present 
work with parameter set 8 are compared with the static nonrelativistic quark model 
(NQM), QCD sum rule (SR) lattice simulations (Latt.) fS^, bag model (Bag) |8l[, 
and Skyrme model (Skyr.) [^2|, The magnetic moments are given in units of the 
nuclear magneton. 



Particle Experiment [H Set 8 NQM SR Latt. Bag Skyr. 



/x(A) 

t(Ap) 
t(S-A) 



2.79 ±10-'^ 


2.82 


2 


82 


3.04 


2.3 


2.78 


1.97 


-1.91 ±10-^ 


-1.66 


-1 


88 


-1.79 


-1.3 


-1.83 


-1.24 


2.42 ± 0.05 


2.48 


2 


70 


2.73 


1.9 


2.65 


2.25 


-1.157 ± 0.025 


-1.09 


-1 


05 


-1.26 


-0.88 


-1.40 


-0.88 


-0.613 ± 0.004 


-0.64 


-0 


60 


-0.50 


-0.41 


-0.60 


-0.59 


-1.250 ± 0.014 


-1.28 


-1 


43 


-1.32 


-0.96 


-1.40 


-1.42 


-0.679 ± 0.031 


-0.78 


-0 


49 


-0.93 


-0.42 


-0.53 


-0.40 


1.261 ± 0.004 


1.25 


1 


67 






1.224 


0.61 


0.718 ± 0.015 


0.760 


1 


00 






0.757 




-0.340 ± 0.017 


-0.238 


-0 


33 






-0.252 




0.250 ± 0.042 


0.190 
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0.167 




1.287 ± 0.158 


1.13 


1 


67 






1.256 





Chapter 6 

Discussion and conclusions 



We have considered a relativistic model for the three-quark core of the baryons. A field- 
theory calculation of matrix elements between bound states is given with the help of the 
quasipotential reduction of the Bethe-Salpeter equation. 

The input to our model fits is the constituent quark masses and the momentum range 
parameters. They are essentially free parameters in the framework of spectroscopic mod- 
els, but they are kept fixed for the entire set of reactions. This is not an easy task because 
the minimal model will not do the job. The physical picture that emerges in our analysis 
is an asymmetric three quark state with a spin-isospin-0 diquark for the hyperons and 
a symmetric wave function for the nucleons, which is Lorentz shaped in the momentum 
distribution. Only for the strangeness-changing weak decay do we need nontrivial form 
factors. With this input we can fit almost every measurement within 6%; many fits are 
even more accurate (see Table |5]^) . 

There is no need for quark form factors in the electromagnetic sector except for fine 
tuning. Wave functions can be described without admixtures of mixed permutation sym- 
metry. In fact, a slight asymmetry in the wave function does correspond to a nonvanishing 
mixed permutation symmetry, but we need a large diquark clustering, as we have seen in 
the previous Chapter. 

The symmetry breaking scheme and the My and of our model can explain all 
present data for the weak beta decay of the hyperons within the Cabibbo theory (see 
Sec. [4.4|) . The value for Vus = 0.225 ± 0.003 that we get has recently been confirmed in a 
different analysis 0]. But a discrepancy with the meson sector still remains. 

We list some features of our model, which have not been looked at before in this 
framework: 

1. Asymmetric wave functions. 

2. Wave functions with admixture of mixed permutation symmetry. 

3. Wave functions with realistic high- momentum features up to more than 30 GeV^. 

4. Comprehensive calculations for both the electromagnetic and the weak sector. 

5. Consistent symmetry breaking scheme for the Cabibbo theory. 
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6. Derivatives of the weak form factors. 

It is interesting to note that the mass parameters in this work are similar to those 
in the meson sector of the same model |T^. This gives us confidence concerning the 



generality of the model. 

We conclude with some remarks on how to improve the model described in this thesis. 

1. In addition to the three valence quarks, one should consider the effects of higher 
Fock states (valence quarks, gluons). 

2. The momentum wave function should be derived from a potential. 

3. Gluon corrections should be calculated. 

4. A different diquark clustering should be investigated. 

Points 1 and 3 have been considered to some extent in our model because diquark 
clustering simulates gluon effects and these also give rise to higher Fock states. Point 2 
should justify the choice of our wave function. 



Appendix A 
Computational methods 



The main programs for this thesis are written in FORM Mathematica FOR- 
TRAN and C [0], together with the NAG library [||]. They run on an Alhant 
FX/80, a Sun SPARC station, a Macintosh Classic, and a Macintosh Ilsi. The final 
version of the programs are written in Mathematica and C to meet two goals: 

1. The software should be easy to modify and to maintain since we want to check many 
different fiavors of the quark model. 

2. Because the formulae are large, there should be a way to check their correctness 
easily. 

The general design philosophy is given in Fig. [A.l| . The input, the description of 
the process, should be small and similar to the physical notation. The properties of the 
Melosh transform and of the various wave functions should be kept in a database, so 
that they have to be typed only once and can be thoroughly checked. From this starting 
point, the symbolic program should produce the large formulae in an automated way and 
splice them into the C program. With this method it is easy to achieve the two above 
mentioned goals, because the input is short and can be given in a physical language. 

The rest of Appendix ^ is organized as follows: In Sec. A.l we give details of the 



symbolic implementation. Numerical questions are treated in Sec. |A.2| , in Sec. |A.tj| we 
show our procedure for the multidimensional fit, and in the last Section some parts of the 
program are listed. 
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Description of process 




C Program"^ 



Fit 



Figure A.l: Flow chart of the basic steps for this thesis. The aim was to keep the input 
for the programs as small and simple as possible to make sure to get correct results. The 
steps in between have been automated to minimize typographical errors. 



A.l Symbolic calculation 



In implementing the wave functions of Eqs. (|2.34|) , (|5.2|) , and (|5 
relations between them. For instance: 

1^) = ~\p) {u d) , 
|E+) = -\p){d^s) , 

\vi) = -bT)(T^i) • 



we use some simple 



(A.l) 



The wave functions in the Mathematica program are written in an obvious notation. 
The flavors up, down, and strange are denoted by their initials u, d, s. The functions in 
Eq. ( |5.4| ) are pi, p2, and p3, respectively. Spin up and spin down are denoted by a and 
b respectively. This part of the Mathematica program reads: 

prot = -(uud (pi (aba-aab)+p2 (baa-aab) ) 
+udu (pi (aab-aba)+p3 (baa-aba) ) 
+duu (p2 (aab-baa) +p3 (aba-baa) ) ) /Sqrt [18] ; 
lam = (sud pl(aba-aab)+usd p2 (baa-aab) 
+sdu pl(aab-aba)+uds p3 (baa-aba) 
+dsu p2 (aab-baa) +dus p3 (aba-baa)) /Sqrt [12] ; 



ab = {aab -> -bba,aba -> -bab,baa -> -abb}; 
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X 



X 



ud = {uud -> -ddu,udu -> -dud,duu -> -udd}; 

us = {uud -> -ssd,udu -> -sds,duu -> -dss}; 

ds = {uud -> -uuSjUdu -> -usu,duu -> -suu}; 

udds = {uud -> ddSjUdu -> dsd,duu -> sdd}; 

usdu = {uud -> ssu,udu -> sus,duu -> uss}; 
Sigma = {uud->(uds+dus)/Sqrt [2] ,udu->(usd+dsu)/Sqrt [2] , 
duu-> (sud+sdu) /Sqrt [2] } ; 

proton [spin_] := If [spin==l/2 ,prot ,prot /. ab] ; 
neutron [spin_] := If [spin==l/2,prot /. ud, 

prot / . ud / . ab] ; 
sigmaP [spin_] := If [spin==l/2 ,prot /. ds, 

prot / . ds / . ab] ; 
sigmaO [spin_] := If [spin==l/2 ,prot /. sigma, 

prot /. sigma /.ab]; 
sigmaM [spin_] := If [spin==l/2 ,prot /. udds, 

prot / . udds / . ab] ; 
;iO[spin_] := If [spin==l/2,prot /. usdu, 

prot / . usdu / . ab] ; 
„iM[spin_] := If [spin==l/2 ,prot /. us, 

prot / . us / . ab] ; 
lambda[spin_] := If [spin==l/2 , lam, lam /.ab]; 

The wave function of the IH*^ I) baryon for example can now be called by xiO [-1/2] . 
The transition matrix elements for the weak decay in Eq. ( |4.7| ) can be defined by a function 
transition [.,.]. 

transition [b2_ ,bl_] := 

Block [{b},b=b2 /. bar; 

Expand [3 b bl] /. flavourWeak /. spin 

] 

This function needs some additional definitions for the spin and flavor part. 

flavourWeak ={udu udd -> l,dud duu -> l,ssd ssu -> 1, 
dsd dsu -> l,sdd sdu -> l,usd usu -> 1, 
sud suu -> l,usu usd -> 1, 
uds udu -> l,dus duu -> l,dds ddu -> 1, 
sds sdu -> l,dss dsu -> l,sus suu -> 1, 
uss usu -> l,sdu sdd -> l,dsu dsd -> 1, 
uud->0 , udu->0 , duu->0 , ddu->0 , dud->0 , udd->0 , 
uds->0 , dus->0 , usd->0 , dsu->0 , sud->0 , sdu->0 , 
uus->0 , usu->0 , suu->0 , dds->0 , dsd->0 , sdd->0 , 
ssu->0 , sus->0 , uss->0 , ssd->0 , sds->0 , dss->0} ; 



spin = {aabs aab->aabaab,aabs aba->aababa,aabs baa-> 
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aabbaa,abas aab->abaaab , abas aba->abaaba, 
abas baa->ababaa,baas aab->baaaab, 
baas aba->baaaba,baas baa->baabaa}- ; 
bar = {aab->aabs , aba->abas ,baa->baas ,pl->pls , 
p2->p2s ,p3->p3s} ; 

We now turn to the spin part of the bracket and write the Melosh transform in 
Eq. ( ^J2^ as: 



Melosh=-[{{a c-qr Ql,-a ql-c Ql}-,{c Qr +a qr,a c-ql Qr}}, 
{■[a d+qr Ql,a ql-d Ql},{d Qr-a qr,a d+ql Qr}}, 
{{b,Ql},{-Qr,b}}} ; 

aab = {Melosh[[l]] .{l,0},Melosh[[2]] .{1,0}, 

Melosh [[3]] .{0,1}}; 
aba = {Melosh [[1]] .{1,0}, Melosh [[2]] .{0,1}, 

Melosh [[3]] .{1,0}}; 
baa = {Melosh[[l]] .{0,l},Melosh[[2]] .{1,0}, 

Melosh [[3]] .{1,0}}; 
abb = {Melosh[[l]] .{l,0},Melosh[[2]] .{0,1}, 

Melosh [[3]] .{0,1}}; 
bab = {Melosh[[l]] .{0,l},Melosh[[2]] .{1,0}, 

Melosh [[3]] .{0,1}}; 
bba = {Melosh[[l]] .{0,l},Melosh[[2]] .{0,1}, 

Melosh [[3]] .{1,0}}; 

For calculating the spin bracket we define the function bracket [.,.], e.g. the transition 
(tit l^^l Tit) can be computed by bracket [aba, aba] . 

o = {{{1,0}, {0,1}}, {{1,0}, {0,1}}, {{1,0}, {0,-1}}}; 
qQrulel={qr->qt2/ql , qr->Qt2/Ql} ; 
qQrule2={ql~n_Integer?Negative-> (qt2/qr) "n, 
Ql"n_Integer?Negative-> (Qt2/Qr) "n} 
qQrule3={ql Qr->qtQt,(ql Qr) "2->qlQr2 , (ql Qr) "3->qlQr3 , 
qr Ql->qtqt,(qr Ql) "2->qlQr2 , (qr Ql) "3->qlQr3} ; 
conjugate={qr->ql , ql->qr , Qr->Q1 , Ql->Qr , a->as , b->bs} ; 
bracket [x_ ,y_] := 

Block [{z},z = X /. conjugate; 

Factor [ReplaceAll [ReplaceAll [ReplaceAll [Expand [ 
Product [z[[i]] .o[[i]] .y[[i]] , {1,1,3}]] ,qQrulel] , 
qQrule2] ,qQrule3]] 

] 

The computation of bracket [aba, aba] takes 40 seconds on a Macintosh Ilsi for = 0: 



In = bracket [aba, aba] //Timing 
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Out = {39.23333333333333333*Second, 
(b*bs - Qt2)*(a"2*as"2*c"2*d~2 + 

2*a~2*c*d*qlQr2 - 4*a*as*c*d*qlQr2 + 
2*as~2*c*d*qlQr2 + a"2*as"2*c~2*qt2 + 
a~2*as"2*d"2*qt2 + a"2*as"2*qt2~2 + 
2*a*as*c"2*d"2*Qt2 + 2*a*as*c~2*qt2*Qt2 - 
2*a~2*c*d*qt2*Qt2 + 4*a*as*c*d*qt2*qt2 - 
2*as"2*c*d*qt2*Qt2 + 2*a*as*d~2*qt2*Qt2 + 
2*a*as*qt2"2*Qt2 + c"2*d"2*Qt2"2 + 
c~2*qt2*Qt2~2 + d"2*qt2*qt2"2 + qt2"2*Qt2~2 
)} 



In order to obtain a convenient form for the numerical implementation in Sec. |A.2| we 
define the function list, which gives a list of the coefficients of the spin parts. 

coef f [a_ ,b_] : =Block [{x} ,x=Coef f icient [a,b] ; 

{Coefficient [x,pl*pls] , Coefficient [x,pl*p2s] , 
Coefficient [x,pl*p3s] , Coefficient [x,p2*pls] , 
Coefficient [x,p2*p2s] , Coefficient [x,p2*p3s] , 
Coefficient [x,p3*pls] , Coefficient [x,p3*p2s] , 
Coefficient [x,p3*p3s] }] 
list [a_] :=Block[{x},x=Expand[N[a]] ; 

Flatten [-[coef f [x,aabaab] ,coeff [x,aababa] ,coef f [x,aabbaa] , 
coef f [x , abaaab] , coef f [x , abaaba] , coef f [x , ababaa] , 
coef f [x , baaaab] , coef f [x , baaaba] , coef f [x , baabaa] }] ] 
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A. 2 Numerical calculation 

The integrals to be handled are normally six dimensional. In some special cases, such as 
for = and for a symmetric wave function, we have only four dimensional integrals to 
do. But it is obvious that a fast and accurate integration routine is crucial for our analysis. 
After testing various library and self- written routines we use the routine dOlf cf from the 
Nag library ||8^. dOlf cf is based on the half subroutine and uses the basic rule 



described by Ref . [0 . The routine attempts to evaluate a multidimensional integral, with 
constant and finite limits, to a specified relative accuracy, using a repeated subdivision of 
the hyper-rectangular region into smaller hyper-rectangles. In each subregion, the integral 
is estimated using a seventh-degree rule, and an error estimate is obtained by comparison 
with a fifth-degree rule, which uses a subset of the same points. The error estimate 
is therefore less time consuming than for the Gauss integration. The fourth differences 
of the integrand along each coordinate axis are evaluated, and the subregion is marked 
for possible future subdivision in half along that coordinate axis which has the largest 
absolute fourth difference. 

For a multidimensional integral it is important to keep the dimensionality of the in- 
tegrals as small as possible. By introducing cylindrical coordinates we can easily reduce 
the six dimensions to five. Defining new variables by 

qi = q cos 9 , Qi = Q cos (j) , 
q2 = qsinO , (^2 = <5 sin , 

q3 = q3 , Q3 = Q3 , 

0<e<27r, O<0<27r, (A.2) 



we can write the difficult parts of the formulae [e.g. ( |5.5|) and ( ^.8|) ] as follows 



q±Q± = qQ cos {6 - (j)) , 
(qiquf = q'Q' cos [2(6 -(P)], (A.3) 
(qiquY = g'Q'cos(^^-0){2cos[2(0-0)]-l} . 

Since the integrand depends only on {6 — (p), one integral is trivial. This is crucial because 
the routine dOlf cf looses much efficiency for more than five integrations. In some cases 
the {6 — 0)-integration can also be done yielding modified Bessel functions of the first 
kind, but the speed of the routine dOlf cf would not increase. 
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A. 3 Fitting procedures 

The procedure to fit our variables to the experimental data is rather involved since the 
number of variables is high for some versions of our model. One way to go is to fit some 
subset of the experiments, e.g. the nucleon sector, and fix the values of the parameters 
so obtained when fitting the rest of the data. The other way is to define a function that 
measures the error of the fit. The Mathematica function for the magnetic moments of 
the octet could be written as: 

f [mu_ , bQs_ , bqs_ , bQx_ , bqx_] : = 
Abs [1 . 42-sP [mu , bQs , bqs] ] + 
Abs [-0 . 16-sM [mu , bQs , bqs] ] + 
Abs [-0 . 61-nil [mu,bQs ,bqs] ] + 
Abs [- 1 . 25-xO [mu , bQx , bqx] ] + 
Abs [0 . 32-xM [mu , bQx , bqx] ] ; 

The calls sP, sM, ml, xO, and xM give the magnetic moments of the S+, S°, A, rP, and 
S~ respectively as a function of the mass of the u quark (mu), and the scale factors /3qE 
(bqs), (bQs), (bqx), and Pqs (bQx). Because of speed considerations we do not 
directly compute the magnetic moments at runtime, but we build a lookup table, which 
we interpolate with the Mathematica function Interpolate, e.g. 

sP=lnterpolation [data_sP , lnterpolation0rder->3] . 

A local minimum of the function f can be found by the MATHEMATICA command 

FindMinimum [f [mu , bQs , bqs , bQx , bqx] , 

{mu,-[.25, .26}, .21, . 33} , {bQs , { . 5 , .52}, .3, .7}, 
{bqs, {.5, .52}, .3, . 7} , {bQx , { . 5 , .52}, .3, .7}, 
{bqx, {.5, .52}, .3, .7}] 

To search for a global minimum we use the graphics capability of Mathematica. 
The command 

DensityPlot [f [.25, .35, .65, bQx, bqx] ,{bQx, .3, .7}, {bqx, .3, .7}] 



produces the plot in Fig. |A.2| . The black areas show the minimum, the white areas the 
maximum of the function f mentioned above. We can easily see in Fig. A. 2 that an 
asymmetric wave function is strongly favored. 
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Figure A. 2: Example of a density plot for the magnetic moments of the baryon octet. 
The black areas show the minimum, the white areas the maximum of the deviation from 
the experimental data. The /5s are given in units of GeV. One can easily see that an 
asymmetric wave function is strongly favored. 




A. 4 Program listing 

In this Section we present some parts of the listing of the C program for the strangeness- 
changing weak semileptonic decay with no momentum transfer {K'^ = 0). The programs 
for 7^ are larger, and the one for electromagnetic properties is shorter. We will not 
present the entire program, but simply comment on some of its parts. The aim was to 
keep the program modular and flexible. 

The FORTRAN routine dOlf cf can be called from within a C program in the fol- 
lowing way: 

dOlf cf _ (fendim , a , b , feminpts , femaxpts , glsn_ , &eps , &acc , felenwrk , 
wrkstr,&suni_glsn,&if ail) ; 

The routine names must be extended by an underscore, and the variables must be called 
by addresses. In the above call we evaluate the form factor gi for the process n 
and store its value in the variable sum_glsn. The function glsn is given below, and the 
parameters are fixed by experience as follows: 

#define LENWRK 100000 

#define N 5 

#define P12 6.283185308 

int ndim=N,ifail=l,niinpts=0,maxpts=(int) (93*LENWRK/7) ,lenwrk=LENWRK; 
static double a [N] ={0 . , -3 . , . ,-3 . ,0 . } ,b [N] ={3 . ,3 . ,3 . ,3 . ,P12} ; 
double eps,acc,wrkstr [LENWRK] , ... 

The succeeding parts of the C program are directly produced by the Mathematica 
programs explained in Sec. |A.1|. The Mathematica command 



FactorTerms [transition [sigmaM [ . 5] , neutron [ . 5] ] ] //list 
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produces the array in the function glsn. 

double glsn_(ndiin,z) 
int *ndim; 
double *z ; 
{ 

double trans [81] ={-0.1666666666666666667, -0.1666666666666666667, 0, 
-0.1666666666666666667, -0.1666666666666666667, 0, 0, 0, 
0, 0.1666666666666666667, 0.1666666666666666667, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0.1666666666666666667, 
0.1666666666666666667, 0, 0, 0, 0, 0.1666666666666666667, 
0, 0, 0.1666666666666666667, 0, 0, 0, 0, 0, 
-0.1666666666666666667, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
-0.1666666666666666667, 0, 0, 0, 0, 0, 0, 
0.1666666666666666667, 0, 0, 0.1666666666666666667, 0, 0, 
0, 0, 0, -0.1666666666666666667, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, -0.1666666666666666667, 0, 0, 0, 0},axial(); 

return (axial (z, trans)) ; 

} 

The function axial, called on the last line, contains the entire implementation of the 
Melosh transform, the kinematics, and the variable transform discussed in the proceeding 
Section. The wave functions are taken from Eq. ( |5.5| ). The spin transitions aabaab, 
aababa, . . . are calculated with Mathematica as shown in Sec. |A.l: 



aabaab=bracket [aab , aab] ; 
aababa=bracket [aab , aba] ; 

Splice ["weak.mc"] 

The last command Splice splices MATHEMATICA output into an external file. The output 
looks as follows: 

double axial (z,t) 
double *z,*t; 
{ 

double q, qt2 , Q , Qt2 , theta , M , M3 , xi , eta , pi , p2 , p3 , pis , p2s , p3s , 

alpha2 , beta2 , alpha2s , beta2s , ml2 , m22 , m32 , m32s , aabaab , aababa , 
aabbaa , abaaab , abaaba , ababaa , baaaab , baaaba , baabaa , res , qtQt , 
qlQr2 , qlQr3 , Ms , Q3 , Q3s , Q_2 , Q_2s , e3 , e3s , e 12 , el2s , a , as , b , bs , c , d , 
Xl,X2,X3,Xls,X2s,X3s,q_2,el,e2,q3; 

q=z[0]; q3=z[l]; Q=z [2] ; Q3=z [3] ; theta=z[4]; 



qt2 = q*q; Qt2 = Q*q; 

q_2 = qt2 + q3 * q3; Q_2 = Qt2 + Q3 * Q3; 

qtQt = q*Q*cos (theta) ; 

qlQr2 = q*q*Q*Q*cos (2 . *theta) ; 
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qlQr3 = q*q*q*Q*Q*Q*cos(theta)*(2.*cos(2.*theta)-l) ; 

el = sqrt(q_2 + ml * ml); e2 = sqrt(q_2 + m2 * m2) ; M3 = el + e2; 

e3 = sqrt(Q_2 + mS * m3) ; el2 = sqrt(C)_2 + M3 * MS) ;M = eS + el2; 

xi = (el + q3)/M3; eta = (el2 + Q3)/M; 

Ms = sqrt((Qt2+(l-eta)*M3*M3+eta*m3s*m3s)/eta/(l-eta)) ; 

qSs = (eta-0.5)*Ms-(M3*M3-m3s*m3s)/2.0/Ms; 

Q_2s = Qt2 + Q3s * Q3s; 

e3s = sqrt(Q_2s + m3s * m3s) ; el2s = sqrt(Q_2s + M3 * M3) ; 

a = M3 + eta * M; as = M3 + eta * Ms; 

b = m3 + (1 - eta) * M; bs = m3s + (1 - eta) * Ms; 

c = ml + xi * M3; 

d = m2 + (1 - xi) * M3; 

alpha2 = alpha*alpha; beta2 = beta*beta; 

alpha2s = alphas*alphas ; beta2s = betas*betas; 

ml2=ml*ml; m22=m2*m2; m32=m3*m3; m32s=m3s*m3s; 

X3= (Qt2/eta/ ( 1-eta) /alpha2+qt2/eta/xi/ ( 1-xi) /beta2+ 

ml2/eta/xi/beta2+m22/eta/ ( 1-xi) /beta2+m32/ (1-eta) /alpha2) ; 

X2=qt2* ( (1-eta) * (1-xi) *alpha2+xi*beta2) /alpha2/beta2 
/eta/xi/ ( 1-xi) /(l-eta+xi*eta)+Qt2*(( 1-xi) 
* (1-eta) *beta2+xi*alpha2)/alpha2/beta2/eta 
/ (1-eta) / ( l-eta+xi*eta) +2 . * (alpha2-beta2) /alpha2 
/beta2/eta/ ( l-eta+xi*eta) *q*Q*cos (theta) + 
ml2/eta/xi/beta2+m22/eta/ (1-xi) /alpha2+m32/ (1-eta) /beta2 ; 

Xl=qt2* ( (1-xi) *beta2+xi* (1-eta) *alpha2) /alpha2/beta2 

/eta/xi/ ( 1-xi) / ( l-xi*eta) +Qt2* ( ( 1-xi) *alpha2+xi* (1-eta) * 
beta2) /alpha2/beta2/eta/ (1-eta) / (l-xi*eta) -2 . * (alpha2-beta2) 
/alpha2/beta2/eta/ ( l-xi*eta) *q*Q*cos (theta) + 
ml2/eta/xi/alpha2+m22/eta/ ( 1-xi) /beta2+m32/ (1-eta) /beta2 ; 
#ifdef GAUSS 

p3=N3*exp(-X3/2.) ; p2=N2*exp(-X2/2 . ) ; pl=Nl*exp(-Xl/2. ) ; 
#else 

p3=N3*pow( (X3+1) , -4) ; p2=N2*pow( (X2+1) , -4) ; pl=Nl*pow( (Xl+1) , -4) ; 

#endif 

X3s= (Qt2/eta/ (1-eta) /alpha2s+qt2/eta/xi/( 1-xi) /beta2s+ 

ml2/eta/xi/beta2s+m22/eta/ (1-xi) /beta2s+m32s/ (1-eta) /alpha2s) ; 

X2s=qt2* ( (1-eta) * (1-xi) *alpha2s+xi*beta2s) /alpha2s/beta2s 
/eta/xi/ (1-xi) / (l-eta+xi*eta) +Qt2* ( (1-xi) 
* (1-eta) *beta2s+xi*alpha2s)/alpha2s/beta2s/eta 
/ (1-eta) / (l-eta+xi*eta) +2 . * (alpha2s-beta2s) /alpha2s 
/beta2s/eta/ (l-eta+xi*eta) *q*Q*cos (theta) + 

ml2/eta/xi/beta2s+m22/eta/ ( 1-xi) /alpha2s+m32s/ (1-eta) /beta2s ; 
Xls=qt2* ( (1-xi) *beta2s+xi* (1-eta) *alpha2s) /alpha2s/beta2s 
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/eta/xi/ ( 1-xi) / (l-xi*eta) +Qt2* ( (1-xi) *alpha2s+xi* ( 1-eta) * 
beta2s) /alpha2s/beta2s/eta/ (1-eta) / ( l-xi*eta) -2 . * (alpha2s-beta2s) 
/alpha2s/beta2s/eta/(l-xi*eta)*q*Q*cos(theta)+ 
ml2/eta/xi/alpha2s+m22/eta/ (1-xi) /beta2s+m32s/ (1-eta) /beta2s ; 
#ifdef GAUSS 

p3=N3*exp(-X3s/2.); p2=N2*exp(-X2s/2 . ) ; pl=Nl*exp(-Xls/2 . ) ; 

#else 

p3=N3*pow((X3s+l) ,-4) ; p2=N2*pow( (X2s+1) ,-4) ; pl=Nl*pow( (Xls+1) ,-4) ; 
#endif 

aabaab=(Qt2 - b*bs) * (Qt2*Qt2*c*c*d*d + 2*Qt2*a*as*c*c*d*d + 
a*a*as*as*c*c*d*d - 2*a*a*c*d*qlQr2 + 4*a*as*c*d*qlQr2 - 
2*as*as*c*d*qlQr2 + Qt2*Qt2*c*c*qt2 + 2*Qt2*a*as*c*c*qt2 + 
a*a*as*as*c*c*qt2 + 2*Qt2*a*a*c*d*qt2 - 

4*Qt2*a*as*c*d*qt2 + 2*Qt2*as*as*c*d*qt2 + Qt2*Qt2*d*d*qt2 + 
2*Qt2*a*as*d*d*qt2 + a*a*as*as*d*d*qt2 + Qt2*Qt2*qt2*qt2 + 
2*Qt2*a*as*qt2*qt2 + a*a*as*as*qt2*qt2) ; 
aababa=(a - as)*(b + bs)*(Qt2*Qt2*c*c*d*d + Qt2*a*as*c*c*d*d + 

28 lines omitted 

These formulae are given in Eqs. (B.41) - (B.57) 

baabaa = -((Qt2 - b*bs)*(Qt2*Qt2*c*c*d*d + 2*Qt2*a*as*c*c*d*d + 
a*a*as*as*c*c*d*d + 2*a*a*c*d*qlQr2 - 4*a*as*c*d*qlQr2 + 
2*as*as*c*d*qlQr2 + Qt2*Qt2*c*c*qt2 + 2*Qt2*a*as*c*c*qt2 + 
a*a*as*as*c*c*qt2 - 2*Qt2*a*a*c*d*qt2 + 
4*Qt2*a*as*c*d*qt2 - 2*Qt2*as*as*c*d*qt2 + 
Qt2*Qt2*d*d*qt2 + 2*Qt2*a*as*d*d*qt2 + a*a*as*as*d*d*qt2 + 
Qt2*Qt2*qt2*qt2 + 2*Qt2*a*as*qt2*qt2 + a*a*as*as*qt2*qt2) ) ; 

res = t [0] *aabaab*pl*pls + t [1] *aabaab*pl*p2s + t [2] *aabaab*pl*p3s + 
t [3] *aabaab*p2*pls + t [4] *aabaab*p2*p2s + t [5] *aabaab*p2*p3s + 
t [6] *aabaab*p3*pls + t [7] *aabaab*p3*p2s + t [8] *aabaab*p3*p3s + 

22 lines omitted 

t [75] *baabaa*p2*pls + t [76] *baabaa*p2*p2s + t [77] *baabaa*p2*p3s + 
t [78] *baabaa*p3*pls + t [79] *baabaa*p3*p2s + t [80] *baabaa*p3*p3s ; 

res /= (a*a+C!t2)*(d*d+qt2)*(as*as+C)t2)*(c*c+qt2)* 
sqrt (b*b+Qt2) *sqrt (bs*bs+Qt2) ; 

res *= sqrt(e3s*el2s*M/e3/el2/Ms) ; 

return (q*Q*res) ; 



The C preprocessor directives #ifdef GAUSS, #else, #endif give us the possibiUty of 
choosing either the Gauss shaped wave function or the Lorentz shaped one by simply 
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typing the line 

#define GAUSS 

or by commenting it out. 

The weak vector form factor is calculated in an analogous way; the line 

o = {{{1,0},{0,1}},-C{1,0},{0,1», {{1,0}, {0,-1}}}; 

on page has just to be replaced by 



o = {{{1,0}, {0,1}}, {{1,0}, {0,1}}, {{1,0}, {0,1}}}; 

The listings for the calculation for ^ are longer, so we omit them here. 
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